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ABSTRACT 

Aims. Spectroscopic observations of Type la supernovae obtained at the New Technology Telescope (NTT) and the 
Nordic Optical Telescope (NOT), in conjunction with the SDSS-II Supernova Survey, are analysed. We use spectral 
indicators measured up to a month after the lightcurve peak luminosity to characterise the supernova properties, and 
examine these for potential correlations with host galaxy type, lightcurve shape, colour excess, and redshift. 
Methods. Our analysis is based on 89 Type la supernovae at a redshift interval z = 0.05 — 0.3, for which multiband 
SDSS photometry is available. A lower-z spectroscopy reference sample was used for comparisons over cosmic time. 
We present measurements of time series of pseudo equivalent widths and line velocities of the main spectral features in 
Type la supernovae. 

Results. Supernovae with shallower features are found predominantly among the intrinsically brighter slow declining 
supernovae. We detect the strongest correlation between lightcurve stretch and the Si n A4000 absorption feature, 
which also correlates with the estimated mass and star formation rate of the host galaxy. We also report a tentative 
correlation between colour excess and spectral properties. If confirmed, this would suggest that moderate reddening of 
Type la supernovae is dominated by effects in the explosion or its immediate environment, as opposed to extinction by 
interstellar dust. 

Key words, methods: data analysis - techniques: spectroscopic - supernova: general - cosmology: observations 

1. Introduction I Kowalski et~aT1 120081: Ifficken et all 120091: lAmanullah et all 

l2010l ). This is emphasised in the first-year SDSS-II cos- 



Cosmological distance measurements based on Type la su- mo i ogv resu i ts Messier et al.l l2009a: Sol lerman et al.1120091: 
pernovae (SNe > la) led to the discovery of the accelerated | Lam p e j tl et AlpX Two of the major (known) sources 



» , tt . , , , itt: n Ltjaiiiouibi cl ai. juiua . i wu ui liic major iaiiuwiii suuices 

expa nsion of the Universe a bout a decade ago (|Riess et al. | of tematic uncert ainties when using SNe la to measure 

USaa J Perlmutter et al. || l999| ), which requires that darken- cosmol ical distances are the corrections for the colour- 

ergy _ exist. This mysterious, hypothetical energy is one of bri htncss rclation andjtpossible drift with redshift of the 

the biggest puzzles m contemporary cosmology and funda- gN brightness ( e . g . lNordin et afl200«h . These shortcomings 

mental physics. ^ _ are related to our lack of a detailed understanding of the 

With the ever increasing statistical precision on the den- underl in ph ics preC eding and during the explosion, in- 

sity and equation-of-state parameter of dark energy, we are duding thg progenitor system and both the c i rcumste ll ar 



now reaching a point where systematic uncertainties are the and intcrstellar environments 
limiting factors (jAstier et al . 2006; Wood-Vasev et al.l2007t 



* -j-j n , ,. „ , , , t. Optical spectroscopy provides an excellent testbed for 

Based on observations collected at the European . 1 . . „ T ,r~ . , . 

Organisation for Astronomical Research in the Southern understanding SNe la: differences m explosion properties 

Hemisphere, Chile, in the ESO programmes 077.A-0437, 078.A- wlU llkel y modify spectral features. Although detailed 3D 

0325, 079.A-0715 and 080.A-0024. Also based on observations modelling may be needed to extract physical information 

with the Nordic Optical Telescope acquired in the programmes (or parameters) from observations, empirical techniques 

with proposal numbers 34-004, 35-023 and 36-010. could provide important hints for future modelling. For 
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example, comparisons between spectral properties for SNe 
with different host-galaxy type, lightcurve parameters and 
redshifts may be used to infer trends that affect their use 
as distance indicators. 

There are indicat ions of a population drift with redshift 
(jSullivan et al.l 120061 ) . which could be detectable as spec- 
tral evolution if the average (composite) spectra in differ- 
ent redshift bins are compared. Redshift dep endencies in 
compo site spectra have al s o bee n suggested by iFolev et all 
([2008D and lSullivan et al.l (|2009t) . While the metal content 
of the Universe increases with time, it is unclear to what de- 
gree increased metallicity actually propagates through the 
progenitor and explosion mechanisms to the outer ejecta 
of the supernova. If the element distribution of the ejecta 
is affected, this will likely change t he observed spectrum 
(jLentz et al.ll2000HSauer et al.ll2008t) . 

An important question is whether the available low-z 
SN data sets correctly sample the demographics of SNe la 
at cosmological distances, or if there are subtypes of SNe 
not (yet) observed in the smaller local samples, but present 
at higher redshifts. Either case could yield an evolution of 
the average spectrum. Our data set at intermediate red- 
shifts provides a useful sample to close the "gap" in studies 
currently available in the literature. We also study indi- 
vidual spectra and are thus potentially able to disentangle 
shifting population demographics from new subtypes. 

It is also important to explore if there is any relation 
between SN la spectral features and broad-band colours. 



to match the Milky- Way dust extinction law dRiess et al. 


1996 


: Tripp! 11998 


iKrisciunas et al. 


20001: Altavilla et al. 


2004 


: Reindl et al. 


12005: Astier et al. 


2006HGuv et al.ll2007t 


Nobili & Goobarl 


2008). The explanation may be con- 



nected to e.g., inter actions in the circumstellar environ- 
ment ([Goobarl [20081) or intrinsic SN properties, in which 
case a correlation of properties of spectral features and 
broad-band colours could be expected. Thus, comparisons 
of spectra of SNe with different colours may help in the 
understanding of the intrinsic colour dispersion in SNe la 
and allow us to disentangle the various components enter- 
ing the colour-brightness relation and its possible evolution 
with redshift, critical for precision cosmology. 

Relations between lightcurve parameters and 
host galaxy pr op erties h ave been pres en ted r ecently 



Kellv et ap (120101) : ISullivan et all 
iLampeitl et al.l (|2010bf ). These empirical 



f|2010h and 
findings call 

for further scrutiny; one way to do so is through compar- 
isons with properties of SN la spectra. 

Although a number of spectral feature comparisons 



Hook et al.l 20051: Benetti et al.l 20051 


Branch et alJ 20061 


Blondin et al.l 12006 


: iGaravini et al. 


2007at Folev et al. 


20081: iBronder et al. 


2008: lEllis et al. 20081: ISullivan et al. 


2009: [Wane et al. 2009: Bailev et al. 20091). we are still far 



from a complete understanding of the observed variation of 
SN la spectra. Detailed spectral studies are needed in order 
to limit the possible differences between low and high red- 
shift objects, a basic requirement for the use of SNe la as 
distance indicators, and could potentially be used to fur- 
ther sharpen the standarizablc candle through secondary 
brightness indicators. 

During 2006-2007, 169 spectra of SNe la were obtained 
at the New Technology Telescope (NTT) and the Nordic 
Optical Telescope (NOT) in a program designed for spec- 



tral identification of objects detected by the S loan Digital 
Sky Survey II (SDSS-II) Supernova Surve y (jGunn et al.l 
[19981: lYork et all [20001 iFrieman et al|[2008t ). The SDSS-II 
Survey operated as a three- year survey (2005-2007), aiming 
at finding a large number of intermediate-redshift SNe la, 
to be used to estimate cosmological parameters. The search 
algorithm and the procedure for the spec t rosco pic observa- 
tions have been described in ISako et ail (|2008l ) . The first- 
ye ar photometry a nd sp ectro scopy have b een p resented 
in iHoltzman et all ([20081 ) and IZheng et al.l ([20081) . respec- 
tively. 

The NTT/NOT SDSS spectra provide a key opportu- 
nity to study SN properties. First, the SN population is 
drawn from an interesting redshift range, where evolution 
could be expected, yet the SNe are close enough to yield a 
reasonably high S/N. Secondly, this data set is large enough 
to allow statistical tests . This data set is described in detail 
in lOstman et al.l (|2010[) . 

We present quantitative measurements of 89 SN la spec- 
tra from the NTT/NOT samples with good lightcurves and 
low to moderate host-galaxy contamination and compare 
these with samples of nearby SNe la. We focus on poten- 
tial evidence of evolution, but also study correlations with 
lightcurve properties, such as stretch and colour, as well as 
host galaxy properties like stellar mass and star formation 
rate. 

The two main sources of systematic uncertainties in 
spectral studies of SNe la are noise degradation and host 
galaxy contamination. These effects complicate the search 
for potential spectral evolution, and would cause systematic 
errors since they will affect nearby and distant SNe differ- 
ently An unknown systematic bias could be misinterpreted 
as a sign of evolution, or even obscure a real effect. 

The analysis of spectral indicators presented here con- 
sists of several steps: 

— We first compare NTT/NOT spectra with nearby data 
(Sect. [4]). Deviating SDSS SNe, possible signs of evolu- 
tion are collected into a deficit subsample (named so 
since measurements are smaller than average). 

— We then change focus and combine all data in order 
to search for correlations with global SN parameters 
(Sect.©. We use LC parameters (SALT and MLCS2k2) 
and host galaxy properties and search for the epoch 
ranges with the most significant correlations. 

— We finally try to understand the origin of both the 
deficit subsample and the major correlations with global 
parameters (Sect.©. This is done through (i) studying 
the spectral region around 4000-4500 A (rest frame) , (ii) 
a comparison of deficit SNe with normal SNe and (iii) 
a short discussion of host galaxy correlations. 

Through each step of this analysis, focus has been put 
on minimising/ mapping any sort of systematic error and/or 
observer bias. The observations, data reduction and host 
galaxy subtraction methods of the NTT/NOT spectra are 
presented in lOstman et al.l (|2010f) . while comparisons of 
indicator measurements will be presented here. Extensive 
Monte Carlo simulations were run in order to estimate er- 
rors caused by host galaxy subtraction or varying noise lev- 
els. 

The full organisation of this paper is as follows: In 
Sect. [2] we introduce spectral indicators and in Sect. [3] we 
present the data sets used. In Sect. @] the indicator mea- 
surements, as a function of spectral epoch, are displayed for 
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Table 1. Feature boundaries (pEW) 



4000 5000 

Restframe Wavelength (A) 



6000 



Fig. 1. The feature regions used for the measurements of 
spectral indicators are shown for two template spectra at 
day —5 and 15 re lative to lightcurve peak (templates from 
iHsiao et al.ll2007t) . The pseudo continuum (dashed orange) 
and line minimum (marked with blue triangle) are shown 
for two features (feature 4 and 7, respectively). 

both the NTT/NOT and the reference SNe. Scct.[3]contains 
studies of the correlations between spectral indicators and 
lightcurve parameters as well as with host galaxy proper- 
ties. Sect.H]contains a discussion (focused on evolution and 
the Si II A4000 feature) and finally conclusions are given in 
Sect. [7] We present further details of our Monte Carlo tests 
in appendices [A] and [B] All SNe included in this study can 
be found listed in Appendix [Cj 



2. Spectral indicators 

Elements in the SN ejecta will absorb photons originally 
emitted by radioactive material in the inner layers, thus 
causing the typical pattern of "features" visible in SN la 
spectra. In this analysis of spectral features we co ncentrate 
on sev en r egions corresponding to these studied bv lFolateTiil 
(|2004t) and lGaravini et all (j2007al) . In Figured] the features 
are displayed for a typical spectrum at two different ages. 

Each feature is labelled after the ion that normally dom- 
inates the absorption in this region (see Table [1| , but since 
most absorption lines are blends of several lines and should 
not be directly identified with physical properties, these re- 
gions will simply be identified as features 1 to 7. 

With the term "spectral indicator" we refer to a mea- 
surement of a spectral feature of a SN la spectrum. Spectral 
indicators are always measured on rest-frame spectra, and 
all spectra presented here have been de-redshifted. In this 
paper we will discuss two indicators, pseudo Equivalent 
Widths (pEWs) and velocities. Feature 5, shaped like a 
V, has two minima; we use the redder of these as velocity 
indicator. 



2.1. Pseudo equivalent widths 

For astrophysical objects with a well-defined continuum, 
the equivalent width of an absorption feature can be mea- 
sured easily. If, furthermore, the density structure is known 
and the absorption is caused by a single ion, this infor- 
mation can be used to deduce elemental abundances. For 



Feature 


Dominating 


Lower region 


Upper region 




line 


(centre A) 


(centre A) 


fl 


Ca II H&K 


3450 - 3800 


3800 - 4100 


£2 


Si ii A4000 


3800 - 3950 


4000 - 4200 


f3 


Mg ii A4300 


3850 - 4250 


4300 - 4700 


f4 


Fe ii A4800 


4300 - 4700 


4950 - 5600 


f5 


S ii W 


5050 - 5300 


5500 - 5750 


ffi 


Si ii A5800 


5400 - 5700 


5800 - 6000 


f7 


Si ii A6150 


5800 - 6100 


6200 - 6600 



SNe la, the spectra are dominated by wide absorption 
features caused by mixed multiple absorption lines. The 
continuum can thus not be read directly from the ob- 
served data, and the physical interpretation of equivalent 
widths becomes non-trivial. Nevertheless, we can measure 
equivalent widths if an unambiguous (pseudo) contin uum 
can be defined We do this following iFolatellil j2004) and 
iGaravini et al . (2007a). A lower and upper limit is found 
at the peak-flux wavelength within lower and upper wave- 
length regions. These regions, for the features used here, 
are given in Table [TJ The pseudo-continuum is defined as 
the straight line between the flux of these lower and up- 
per limits, with the choice of peaks optimised so that the 
pseudo continuum is maximised while not intersecting the 
spectrum. 

With this pseudo-continuum a (pseudo) Equivalent 
Width (pEW) can be calculated (using the standard equiv- 
alent width formula): 



pEW = 



1 - 



A A,; 



(1) 



where / is the observed flux and f c is the pseudo- 
continuum. The sum is taken over all wavelength bins con- 
tained between the lower and upper limit. 

The pEW definition has the advantage of not having 
to fit any function to the data (most features are clearly 
non-Gaussian) as well as being insensitive to multiplicative 
differences between spectra (assuming the multiplied factor 
does not change drastically over the range of the feature). 
However, pEWs are not insensitive to additive flux differ- 
ences. 

As for equivalent widths, the statistical error is given by 



2 

"pEW 




(2) 



It consists of two parts, the first is obtained from the error 
spectrum, cry, while the second propagates the uncertainty 
from the choice of pseudo-continuum, cr/ e . 

To avoid subjectivity in the pEW measurements, all 
steps were automated, e.g. the level of filtering was deter- 
mined through lookup tables (see below and Appendix B) 
and boundaries for the pEW are determined using com- 
puter algorithms. The automated code w as validated by 
measur ing the pEWs on the same data as IGaravini et all 
(|2007ah . The same indicator trends were obtained when us- 
ing the same input spectraQ 



1 As a further check on these algorithms the automatic mea- 
surements were manually revised and the outcome compared 
with our basic results. No major deviations were seen. 
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In addition to the statistical error, there are several 
sources of systematic uncertainties. Host-galaxy contami- 
nation can both change the shape of the feature and in- 
duce an additive flux change. Differential slit loss effects, 
which although being multiplicative can have large effects 
on the host-galaxy subtraction process, are included in con- 
tamination errors. These systematic pEW uncertainties are 
discussed below. 



Noise-filtering uncertainties Filtering or smoothing is nec- 
essary in order to identify the end points of spectral fea- 
tures. However, the optimal filter parameters will change 
with Signal-to- Noise (S/N) ratio; noisy data need stronger 
smoothening to reduce the impact of random fluctuations 
in choosing the feature endpoints, while the same filter 
strength will dilute information in high S/N data. Over 
or under filtering can introduce a measurement bias. This 
is a particular concern when making redshift comparisons, 
since distant SNe have (in general) lower S/N than nearby 
SNe. 

Extensive Monte Carlo (MC) simulations were carried 
out to make an optimal choice between filtering methods 
and their parameters. The simulations are described in de- 
tail in Appendix [B] We conclude that a standard boxcar 
smoothing performs well, provided that the boxcar width 
is modified depending on the noise level of the spectrum 
and what feature is being studied. The simulations were 
used to find the optimal boxcar width for each S/N. For 
each feature and S/N level we also obtain an uncertainty 
from the simulations which is added in quadrature to the 
systematic pEW error. 





> A- - 




--i-- 




+ 
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E(B-V) (mag) 

Fig. 2. Percent change in pEW for features 2 to 7 as a 
function of E(B— V). The measurements are based on the 
iHsiao et"aL (l2007f) SN temp late at peak brightness and as- 
suming ICardelli et ahl ()1989t ) type dust with R v = 1.7. 

Table 2. Feature minima wavelengths 



Feature 


Dominating line 


Rest wavelength (A) 




(~ A observed) 


fl 


Ca II H&K 


3945.12 


£2 


Si ii (A4000) 


4129.73 


f3 


Mg ii (A4300) 


4481.20 


f4 


Fe II (A4800) 


5083.42 


£5 


S ii W 


5536.24 


ffi 


Si ii (A5800) 


6007.70 


f7 


Si ii (A6150) 


6355.21 



Host-galaxy contamination uncertainties The host galaxy 
contamination is the single largest source of systematic er- 
rors for pEW studies. Unsubtracted host light can both 
cause an error through a flux offset and, if the underly- 
ing galaxy is changing with wavelength, through shifting 
pEW feature boundaries. Monte Carlo simulations were 
performed to estimate the uncertainties due to the host- 
galaxy subtraction. When host galaxy contamination could 
be estimated using photometry, both the error and a possi- 
ble bias is retrieved as a function of contamination (this is 
the case for all NTT/NOT spectra as will be discussed be- 
low). For spectra where no host-galaxy is subtracted (refer- 
ence spectra and low contamination SDSS spectra) we add 
the uncertainty expected for an uncorrected 10% galaxy 
contamination. See Appendix [X] for a more detailed de- 
scription of these simulations. 

Reddening Pseudo-equivalent widths are also affected by 
un-corrected reddening by host-galaxy dust. In Figure [2] 
we show how pEWs, as measured on a template spectrum, 
change as dust-like extinction is applied to the template. All 
features gradually decrease with colour excess, E(B — V). 
Changes are smaller than 10% for E(B — V) < 0.3 mag. As 
expected, wider features change more with extinction than 
narrower ones. 

2.2. Line velocities 

The position of absorption and emission features can also 
be used to probe SN properties. We study the wavelength 



minima of the features defined above. As reference minima 
we use the rest wavelength of the ion that each feature was 
named after, these are given in Table [2] (3rd column). 

These are converted into velocities through the rclativis- 
tic Doppler formula, 

_ (A m /A ) 2 - 1 

V a bs — C , [J) 

{^ m /Xor + 1 

where Ao is the laboratory wavelength of the ions creating 
the feature and X m is the measured wavelength in the rest- 
frame of the host galaxy (SNe without host galaxy redshifts 
are thus excluded from velocity studies). 

As for pEWs, it is difficult to give a direct physical in- 
terpretation of line velocities for SNe la since most spec- 
tral features consist of blends of ions. Also, different ions 
dominate features at different epochs and thus shift the 
minimum position. Nonetheless we use the same reference 
wavelength, with the consequence that measurements are 
not guaranteed to be the velocity of an ion, but are rather 
a general measurement that can be compared between dif- 
ferent SNe. In practice we only study each feature during 
epochs where no drastic changes happen to the minima 
shape. Some of the features lack clear minima, like feature 
4, and are thus harder to evaluate. Si II A6150 (f7), typical of 
SNe la, most often yield unambiguous measurements and 
is thus usually the best estimate of the "true" expansion 
velocity. 

All spectra are binned using bin widths constant in ve- 
locity, cAA/A = 2000 km s" 1 . Binned error spectra were 
calculated through weighted error averages in each bin. 
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The flux minima for the different features were obtained 
in the binned spectra, i.e. without interpolation to a finer 
wavelength grid and without fitting a general shape to the 
region. Although a sub-bin fitting procedure would, in prin- 
ciple, give a finer determination of the minimum, this is not 
viable with the noise level in our sample. 

Velocity error estimates Velocity measurements are less sen- 
sitive to most major systematic uncertainties, but noisy or 
heavily contaminated spectra can still have substantial un- 
certainties. These were studied using similar Monte Carlo 
simulations as for pEWs, see Appendix [Al and IB"1 For noisy 
or contaminated data, the dispersion is non-negligible, but 
leads to no major bias, provided host-galaxy subtractions 
are performed on host contaminated spectra. 

Systematic uncertainties obtained from the simulations 
are added to each velocity measurement depending on the 
S/N and the estimated contamination. These are added in 
quadrature with a 200 km s _1 peculiar velocity error. 

3. Data sets 

Our sample consists of SNe la observed with the NTT 
and the NOT as a part of the SDSS-II Supernova Survey 
(jOstman et al.l[2010t ) . The set consists of 169 SN la spectra 
of 141 different objects. SDSS SNe are labeled after their 
SDSS Supernova ID number, see lOstman et al.l (|2010F) for 
the respective IAU names. 

Observations made at the NTT were performed using 
the ESO Multi-Mode Instrument (EMMI) and have a wave- 
length coverage from 3800 to 9200 A, a wavelength disper- 
sion of 1.74 A per pixel, and a spatial resolution of 0'.'166 per 
pixel before binning. A binning of 2x2 was used. The NOT 
spectra were obtained using the Andalucia Faint Object 
Spectrograph and Camera (ALFOSC) with grism 4. NOT 
spectra have a wavelength range from 3200 to 9100 A, a 
wavelength dispersion of 3.0 A per pixel, and a spatial res- 
olution of 0'.'19 per pixel. See lOstman et al.l (|2010h for de- 
tailed information. 

Lightcurve properties such as stretch, colour and max- 
imum absolute magnitude, as well as the spect ral epochs, 
were obtained with the SALT lightcurve fitter (|Guv et al.l 
l2005h . The spectral epoch is defined with respect to the 
peak of the i?-band lightcurve. 

After applying lightcurve quality cuts, requiring photo- 
metric observations both prior and post maximum bright- 
ness, we are left with 127 spectra. Out of these, 116 spectra 
have both good host-galaxy subtraction and arc of sufficient 
quality for spectral features to be identified. Finally, wc 
apply a host-galaxy contamination cut of < 60% in the g- 
band, motivated by Monte Carlo simulations, which leaves 
us with 89 spectraH A list of all NTT/NOT spectra used 
in this analysis is given in Table IC.21 

The SDSS NTT/NOT spectra are compared to a 
low-redshift reference SN sample which consists of three 
subsets, data from the Harvard- Smithsonian Center for 
Astrophysics (CfA), the Supernova Cosmology Project 
(SCP99) and the Online Supernova Spectrum Archive 



2 This limit is somewhat arbitrary: many subtractions of 
higher contaminated spectra succeed, but the risk of a subtrac- 
tion significantly failing increases above 60% host contamina- 
tion. See Appendix El 



(SUSPECT). Since the NTT/NOT spectra cover the spec- 
tral epochs between -9 days and +20, we have only studied 
reference spectra up to epoch 30. 

T he CfA sample c onsist s of 162 spectra of 19 SNe la 
from iMatheson et al.l (|2008t) . The SCP99 data set contains 
79 spectra of 16 SNe observed by the Super nova Cosmology 
Project in 1999 that were studied by iGaravini et al.l 
(|2007al ). The SUSPECT data set collects publicly available 
SN spectr£0, wc use 421 spectra of 40 Type la Sne. A list 
of all spectra in our reference sample is given in Table IC.ll 
The table also contains the source of the lightcurve pa- 
rameters as well as the original spectroscopic reference for 
SUSPECT spectra. These lightcurve parameters are lack- 
ing for some SNe and these objects are thus excluded from 
analysis where such information is required. 

We present in Figure [3] the distribution of epoch, red- 
shift, SALT stretch and SALT colour for the NTT/NOT 
sample together with the reference sample used in this pa- 
per. 

The NTT/NOT sample has significantly larger rcdshifts 
than the comparison sample, a median value of zsdss = 
0.17 compared to z re f =0.01. The spectral epoch distribu- 
tion is also somewhat different: the NTT/NOT spectra are 
more centred around the lightcurve peak, while the com- 
parison sample includes earlier and later epochs. This is a 
natural effect arising from the differing magnitude limits of 
the SN searches. The distribution of SALT stretch is fairly 
constant between the samples with a Kolmogorov-Smirnov 
(KS) probability of 5%, thus showing that we can not reject 
the assumption that they belong to the same distribution. 
The median value for the NTT/NOT sample and the ref- 
erence sample are 0.96 and 0.94, respectively. However, the 
SALT-c (colour) distributions are significantly different as 
can been see by a visual inspection, with a tail of red colour 
SN in the reference sample. This is not surprising since lo- 
cal SNe la can be detected even with a few magnitudes of 
extinction. The median value for the NTT/NOT sample 
and the reference sample are 0.05 and 0.10, respectively. 

While the SALT lightcurve fit output is used to make 
the initial anal ysis, we have also obtained MLCS2k2 
(jJlia et all [20071) lightcurve fits. MLCS output consist of 
a lightcurve shape-dependent parameter A and the V- 
band extinction, Ay. The Ay paramcterisation assumes 
that any reddening not corrected for by the A parame- 
ter can be described by a Milky Way-like extinction law. 
These fits are used to validate results found using SALT 
and study lightcurve fit dependant effects. We use the 
MLCS lightc u rve fi ts of all reference SNe contained in the 
iHicken et all (|2009D data set together with MLCS fits of 
the NTT /NOT SPS S SNe using the SNANA fit package 
(jKessler et al.ll2009bl) and emp l oying t he same quality cuts 
and settings as in Kessler et al.l (|2009aj) . The only exception 
to this procedure was th at we used Ry = 1.7 , this change 
was made to comply with lHicken et al" (|2009| ) but has very 
small effect on our analysis. 

In lOstman et al.l (|2010l ) three potential peculiar SNe 
la from the NTT/NOT data set are presented: two of 
SN1991T-type and one SN2002cx like. In this paper we 
can not confirm any addi ti onal 'SN1991T', 'SN1991bg' or 
'SN2002cx' SNe. iLi et al.l (|2010f ) find, in their luminosity 
limited sample, that 77% of all Type la SNe arc normal, 
18% SN1991T-like, 4% SN1991bg-like and 1% SN2002cx- 



http : / /bmf ord . nhn . ou . edu/ suspect/ . 
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Fig. 3. Distribution of epoch, redshift, SALT stretch factor and SALT colour for the SN spectra used in our study. The 
epoch is defined here as the number of days in rest frame from B band maximum brightness. The white histogram is 
for the reference sample while the striped histogram is for the NTT/NOT spectra used here. Legends show mean and 
Gaussian la levels for the subsets. 



like. While 'SN1991bg' and 'SN2002cx' cither would have 
escaped detection entirely, due to their low luminosity, or 
have been singled out based on lightcurve properties, it is 
likely that a fraction of the SNe used here are really of the 
'SN1991T'-type. That these are not identified can be ex- 
plained by two effects: (i) The S/N is often good enough to 
classify a SN as la, but not high enough for a strict sub- 
classification ()Ostman et al.|[2010h (ii) 'SN1991T' like SNe 
are often mistaken for n ormal if no early spectra exist (an 
'age bias'- ILTet al.ll2010h . We will later, in Section|6l show 
that a fraction of ~ 20% of the SDSS SNe have shallower 
features and would represent 'SN1991T' like SNe at later 
epochs well. 

The NTT/NOT spectra have well documented uncer- 
tainties (jOstman et al.l[2010T ) . However, most of the spectra 
in the reference sample lack error estimates. For such cases, 
a constant flux error of 5% of the average spectral flux was 
used to compute the uncertainties in spectral indicators. 

3.1. Host galaxy subtraction 

A large fraction of the SDSS NTT/NOT SNe have signifi- 
cant host galaxy contamination which could affect spectral 
indicators. In particular, comparisons with virtually host- 



free spectra for local reference SNe could lead to systematic 
differences that may be confused with evolution with red- 
shift. Great care was taken to remove the host galaxy light 
efficiently, given the data available, as well as understand- 
ing remaining errors. 

A large fraction of the NTT/NOT SNe were not ob- 
served in paralactic angle, and are thus affected by differ- 
ential slit losses. In lOstman et al.l (|2010h we calculate slit 
loss for all spectra, as functions of estimated seeing condi- 
tions and centring wavelength. Because of the uncertain- 
ties associated, mainly regarding centring, we choose not 
to apply calculated slit losses directly to the spectra, but 
rather to incorporate corrections into the host subtraction 
pipeline, where we also account for reddening. Considering 
the typical centring wavelength (6500 A) and the limited 
wavelength range studied (4000-6500 A), slit loss and red- 
dening exhibit similar differential attenuation and can be 
fitted together. 

A detailed d escription of the hos t galaxy subtraction 
can be found in lOstman et al.l (|2010[) . In brief, the galaxy 
SED which is subtracted is estimated by minimising the 
difference between the observed spectrum and a combina- 
tion of a SN template and a set of galaxy eigenspectra. To 
the SN template, a second degree polynomial is multiplied 
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to account for reddening (e.g. due to host galaxy dust ex- 
tinction) and differential slit loss effects. The minimisation 
can be described with the formula 

3 

/fit (A) = a s(A) • / SN (A) + Y, a ^(A), (4) 

?:=i 

where /sn is the SN template, <?i the galaxy eigenspec- 
tra, s the second degree polynomial and dj weights which 
are fitted in the subtracti on. The SN temp lates in the fit 
were the Hsiao templates (iHsiao et al. 2007) within epochs 
±5 days from the SN spectral epoch as obtained from the 
lightcurves. Models of the peculiar SNe 1991bg and 1991T 
(jNueent et al.l [20021) within the same epoch interval were 
also included in the fit, but no new clear cases of these sub- 
types were foun d. Three galaxy eigenspectra from SDSS 
(|Yip et all 12001 were used for the galaxy SED in the fit. 
The second degree polynomial s(A) was locked to have s = 1 
at A = 6600 A and s < 1 for all other wavelengths (thus 
only having one degree of freedom). The wavelength of the 
s function peak is chosen to match the wavelength where 
most spectra were centred on the slit. It should be noted 
that the slit loss function is asymmetric around the centring 
wavelength, but the fit during subtraction is only made be- 
tween 4000 and 6000 A and thus the behaviour of s at longer 
wavelengths will not affect the fit. The polynomial is only 
multiplied with the SN SED, and not with the galaxy. This 
was done since galaxies, not being point sources, are signif- 
icantly less affected by slit loss. 

Several modified versions of this host subtraction were 
tried: Other sets of eigencomponent spectra, more eigen- 
componcnt spectra, free polynomial (instead of restricted) 
and slit loss applied to the host galaxy. For individual spec- 
tra one of these alternative methods might achieve better 
fits, but globally they were all either less stable or as good 
but with significantly more parameters to fit. 

Figure 2] shows subtraction samples of SNe with con- 
tamination ranging from very high (77%) to very low but 
including reddening/slit loss. For moderate contamination, 
10 to 60% galaxy light in the g-band, the multiplicative 
host-galaxy subtraction works well. Within these limits the 
best fit galaxy SED is subtracted from the observed spec- 
trum before any measurements are done. An error associ- 
ated with host contamination uncertainties was calculated 
for each indicator measurement. In Figure [5] we show ex- 
amples of how host subtraction affects pEW measurements. 
For spectra with very low contamination, no subtraction is 
done and the error estimated for 10% unsubtracted light 
added. Spectra with very high contamination are excluded 
from the analysis. 

In Appendix IA1 we describe simulations designed to esti- 
mate uncertainties and evolution detection limits. We have 
performed several tests in order to determine the stability 
of our results. These include comparisons with synthetic 
spectra and comparisons between different host galaxy sub- 
traction methods. 

The low-z reference sample spectra were not host sub- 
tracted. These are sufficiently local to allow subtraction 
of most host galaxy contamination during data reduction. 
Visually, they do not appear to contain significant host 
galaxy light. 



4. Results: Comparing the reference and 
NTT/NOT samples 

After host subtraction all spectra are processed through 
the automated indicator measurement pipeline. The error 
bars of the measurements are symmetric geometric sums 
of the statistical uncertainty of the measured indicator and 
the noise-filtering and host-galaxy subtraction systematic 
errors. 

All spectral epochs used are rest frame epochs and are 
thus corrected for time dilation. 

4.1. The reference sample 

The reference set measurements were combined into l-er 
contours to facilitate statistical comparisons. The contour 
is calculated for each day as the weighted mean and un- 
certainty of the indicator (pEW/velocity) for ±3 days. The 
broad epoch interval is used to make a smooth curve (stable 
with respect to outliers). As a justification for the defini- 
tion of the band, the measurements underlying the band 
for pEW for feature 3 are shown in Figure [5] 

We confirm the displacement of the unusal SNe la com- 
pared to the overall trend, as shown in Figure [7] for the 
reference sample. Peculiar SNe are not included in the ref- 
erence sample, but it is in practise impossible to make a 
strict definition regarding which SNe should be considered 
as "normal" . The number of observed spectra per SN also 
varies, thus giving artificially high weight to certain ob- 
jects. We thus do not expect the reference sample to com- 
pletely match the NTT/NOT SNe detected in a rolling SN 
survey. The fraction of SNe of different SN la subtypes are 
different, as well as the distribution of lightcurve color. We 
will return to possible consequences of this when discussing 
evolution in Section [B] 

4.2. Pseudo-equivalent widths 

Figure [5] shows the measured pEW values of features 2 to 7 
for the NTT/NOT spectra as a function of epoch. These are 
compared to the corresponding 1 a contour for the normal 
SNe la in the reference sample. Features are only measured 
for epochs when they can be clearly defined. 

To be able to study the correlations of spectral indica- 
tors with different parameters we want to remove the epoch 
dependence. This is done by fitting a function describing the 
epoch evolution in pEW and then subtracting it from the 
measurementsQ. This epoch independent quantity (ApEW) 
is shown in Figurc[S]for all SNe la NTT/NOT and reference 
spectra (including peculiar types) as a function of redshift 
for epochs around maximum brightness. 

Wc now examine Figure [H] and Figure [5] for significant 
differences between distant and local SNe. The NTT/NOT 
measurements generally match the l-er contour of the refer- 
ence sample within uncertainties. There are, however, some 
regions where the NTT/NOT pEW measurements appear, 
on average, lower than the reference set average. These dif- 
ferences appear most significent for pEW f2 and pEW f4, 



4 Indicators that vary little or gradually with epoch are fitted 
with a linear function. Indicators that show sudden changes (like 
pEW f3) are fitted using a logistic function: f(t) = A/(l + 

g(*i.r-t)/T) _|_ % 
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Fig. 4. Display of sample host galaxy subtractions of SDSS SN16637, 12907, 17886 and 13894 (clockwise from top left, 
restframe epochs -1, 0, -4 and 9). All flux values (y-axis) have been scaled. Wavelength (x-axis) is rest-frame wavelength. 
The grey line shows the unsmoothed original spectrum, while the dashed blue line shows the smoothed version. Dotted 
red line show final subtracted slit loss/reddening corrected spectrum and orange solid line show best fit galaxy. Note that 
because of the slit loss/reddening correction the galaxy and subtracted spectrum do not sum to the raw spectrum. The 
galaxy contamination estimated from photometry (g-band) and redshift have been written in each panel. 
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Fig. 6. pEW values for the reference data for f3 vs epoch. 
The shaded region is the same one sigma contour as shown 
in the left panel of Figure[U The different symbols show the 
measurements that were used to construct the grey region, 
with different symbols denoting the different subsets of the 
reference sample. 



Fig. 7. Peculiar SNe compared with normal. The shaded 
region shows the one sigma contour constructed from the 
normal SNe la in the reference sample and the NTT/NOT 
spectra. The symbols show the measured pEWs for the SN 
1991T-like, SN 1991bg-like and peculiar SNe la in the ref- 
erence sample. 



this is most easily seen in Figure O (top left and mid left, 
shaded region). 



In order to study the origin of these differences we col- 
lect all SNe with cither f2 or pEW measurements be- 
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Fig. 5. Expanded view of the feature 3 region for the same spectra displayed in Figure 0] showing the effects of host 
subtraction and slit loss/reddening correction on pEW. Wavelength (x-axis) is rest- frame wavelength. The blue dashed 
line is the raw spectrum, the red solid line the final subtracted and corrected spectrum. The marked regions show 
where pEW is calculated, the calculated pEWs ( in A) are written in panels for raw spectra ("raw"), subtracted spectra 
( "sub" ) as well as subtracted and corrected spectra ( "slit" ) . These sample spectra show host subtraction change pEWs 
both through additive offsets and changed feature limits. 



low those of normal SNe in the reference sample into a 
pEW-deficit sample (thus ApEW-f2 < -8 A OR ApEW- 
f4 < — 25 A). This subset is examined in detail in Section [51 
where wc for example discuss the effects of different 
lightcurve parameter distributions. The limits in pEW for 
the deficit sample are arbitrary; we do not expect them to 
precisly single out a physically distinct subset of SNe. It 
is rather to be seen as a starting point for a discussion of 
possible differences between local and distant SNe. 



4.3. Comparing the two samples - line velocities 

Figure [TU] shows some of the more stable line velocities for 
the NTT/NOT sample together with the 1-er contour for 
the normal SNe la in the reference sample. Only measure- 
ments of SNe with redshifts measured from galaxy lines are 
included, i.e. the subset of objects with redshifts measured 
from SN features are excluded here. No signs of sample 
differences are detected. 

The epoch evolution of the velocity of f7 Si ii A6f50 
(with epoch) have been extensively st udied ( Bc netti et al.l 
120051: IWang et al.l 120091: iMaeda et al1l2010l) . However, the 
NTT/NOT sample does not contain enough SNe with mul- 
tiple spectra to measure velocity changes. 



4.4. Summary: Comparing the reference and NTT/NOT 
samples 

The samples are generally consistent, possibly deviating in 
a subset of NTT/NOT SNe with pEW measurements below 
these for normal SNe la in the reference sample. These were 
collected in a pEW-deficit sample. 



5. Results: Correlations with SN parameters 

The collected sample (both low- and high-z SNe) was used 
to search for correlations between spectral indicators and 
global properties of Type la SNe. Since many of the features 
evolve with epoch, we study the epoch corrected pEW- and 
velocity-differences, ApEW and Av, as introduced above. 

The correlation is calculated taking into account the 
estimated uncertainties of the indicators. 



5.1. Correlation statistics 

As the basic measure of correlation between measurements 
Ri, Si, (i = l..n), we use Spearman's rank correlation coef- 
ficient, 



r s = 



E t (Rj-RKSj-S) 



(5) 
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Fig. 8. Study of pseudo-equivalent widths for feature 2 to 7 vs epoch. The shaded band shows the one sigma contour 
for the normal SNc la in the reference sample. The points show the measurements on the NTT/NOT spectra, where the 
error bars include both statistical and systematic errors. The differe nt symbols show two ca tegories: spectra identified as 
of the normal SNc la subtype by SNID (SupcrNova IDentification: [Blondin fc TonrvlfeOOTj ) or spectra identified as SNe 
la but of unknown subtype. 



This is similar to the Pearson correlation coefficient, i.e., a 
non-parametric measure of correlation, but relies on ranked 
variables. This method is preferable for variables not follow- 
ing a Gaussian distribution. Spearman's rank correlation is 
also less sensitive to outliers. The output coefficient range 
from —1 to 1, with —1 being perfect negative correlation, 



no correlation and 1 being perfect positive correlation. 
The significance of a correlation r from n elements can be 
estimated roughly using Student's t distribution of dimen- 
sion (n-2) and t = r^/(n - 2)/(l - r 2 ) (jPress et alj|1992f ). 
For example, rank correlation 0.6 corresponds to less than 

1 % chance of being random if n > 18. In our analysis 
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Fig. 9. Comparison of the pseudo-equivalent width measurements around lightcurve peak for features 2 to 7 vs redshift. 
In all figures the average spectroscopic evolution among non-peculiar low-z reference SNe has been subtracted. It is thus 
the pEW- difference (compared to normal low-z SNe) that is plotted on the y-axis. SNe are divided according to subtype 
classification, with the colour scheme following the legend in the upper right plot. In the two upper left panels (f2 and 
f4) the shaded region show how the pEW- deficit sample is defined. 



we mark any correlation r > 0.6 where rt > 15 for further 
study. More realistic confidence analysis should be done 
using permutation tests. This is done when using flexible 
ranges below, where we account for the fact that we probe 
a large number of correlations (and thus expect statistical 
fluctuations to cause some large |r|). We first present the 



basic correlation coefficients for spectra close to lightcurve 
peak. 
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Fig. 10. Comparisons of line velocities (f2, f3, f5, f7) between the reference sample and the higher redshift SDSS spectra. 
The shaded band shows the one sigma contour for the normal SNe la in the reference sample. 



5.2. Correlations with lightcurve parameters 

Correlations around lightcurve peak We have searched for 
correlations with SALT stretch and colour as well as abso- 
lute magnitude, M, for spectra within ±3 days from max- 
imum brightness. The absolute magnitude is corrected for 
stretch and colour and calculated assuming a fiducial cos- 
mology. Table [3] lists the calculated Spearman coefficients. 
The statistical significance in standard deviation from the 
null correlation hypothesis is shown in parenthesis. At least 
50 measurements were used in each correlation estimate. 

Some pseudo-equivalent widths around peak luminos- 
ity do show strong correlation with lightcurve parameters. 
These include f2 and f7 correlating with stretch and f4 
showing a correlation with lightcurve colour. In Figure [11] 
we show these strong correlations, together with the f6 cor- 
relation with stretch. For the latter, all SNe except the low 
S/N NTT/NOT objects show a strong correlation; this fea- 
ture is too small to probe among noisy data. The correlation 
of the depth of this feat ure with lightcurve widt h has been 
reported earlier, see e.g. lHachinger et al.l (|2008h . 

Also noticeable is that most correlations have the same 
direction (sign): SNe with wide lightcurves (large stretch) 
have weaker pEW values and redder SNe (large colour) have 
in general larger equivalent widths. While these colour cor- 
relations are not strong they are consistent and opposite in 



direc tion to what would b e expected from an application of 
pure ICardelli et ail (|1989D type extinction (see Figure [2]). 

Any correlation between spectral indicators and peak 
magnitudes corrected for stretch and colour would be of 
great interest, since this could be used to "sharpen" Type 
la SNe as standard candles. Most spectral indicators corre- 
late only weakly with absolute magnitude (after correction 
for stretch and colour). The most significant correlation is 
for pEW f4; this is, however, a weak correlation of only 
moderate significance (~ 3er). 

We find that none of the velocities are strongly corre- 
lated with stretch, colour or absolute magnitude around 
peak brightness. The kinetic energy, as sampled by line ve- 
locities, thus appear to be independent from the optical 
luminosity. 

Flexible epoch ranges There is no a priori reason to ex- 
pect a fixed epoch range around lightcurve peak to be the 
epoch range where spectral indicators correlate best with 
lightcurve parameters. In Figure [12] we present a sample of 
the bracketed epoch ranges where the most significant cor- 
relations were found through a blind search involving all 
indicators. Before analyzing the results we will describe the 
blind search in detail, as well the MC studies performed in 
order to determine how significant the search output is. 
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Table 3. Correlations at peak brightness: rs, Spearman correlation coefficient (number of standard deviations from 
null). At least 50 measurements used for each entry. 
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Fig. 11. These correlation studies show the measured pEW minus the expected pEW for the corresponding epoch for the 
full sample of normal SNe la vs SALT lightcurve parameter (stretch, colour) for epoch ranges around peak (±3 days) . 
Indicators chosen are the ones showing largest correlation for this epoch range. 



For a given set of epochs, an indicator measurement 
(e.g. pEW f2) and a global property (e.g. stretch), we loop 
through all epoch ranges (containing at least 15 measure- 
ments) and save any correlations with \rs\ > 0.4. A combi- 
nation is thus defined by (spectral indicator, global prop- 
erty, min epoch, max epoch) e.g. ('pEW f2','stretch',-8,0). 
The number of such combinations, for each indicator, range 
from 500-1000 (if using all SNe and an indicator well de- 
fined at all epochs) to 50-100 (if using a subset of SNe and 
an indicator not existing at all epochs). 

For any real correlation found we also expect "neigh- 
bouring" epoch ranges to be correlated. If, for example, 



the epoch range — 8 yields a strong correlation we would 
also expect epoch ranges like 1 — 7 and 2 — 9 to show cor- 
relation. We thus rank correlations between indicators and 
global properties through the number of epoch ranges with 
\rs | > 0.4. When we discuss a correlation in an epoch range, 
this is thus only one in a series of neighbouring correlating 
epoch ranges. 

While this epoch bracketing is necessary in order to find 
the epoch ranges where indicators are sensitive to global pa- 
rameters, this method increases the probability of finding 
random correlations. For any set of three parameters, it will 
always be possible to find some correlation between two of 
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these through restrictions of the third. For any correlating 
parameters we thus have to find the probability of finding 
such corrclation(s) in random data. This is done through 
Monte Carlo simulations where we retain the epoch and in- 
dicator values and randomise the global property. We then 
search for correlations among epoch ranges exactly as for 
real data. This process is repeated 1000 times and the num- 
ber of strongly correlating epoch ranges is saved for each. 
This result can be used to find the probability of finding as 
many strongly correlating epoch ranges by chance. For the 
correlations presented below, we find either zero or one out 
of 1000 iterations to yield as many correlated epoch ranges. 
We thus find that the probability that these correlations are 
completely random is equal to or smaller than 0.001. 

The strongest correlation found was for the pEW for fea- 
ture 2 and stretch when probing epochs right before max- 
imum brightness. Depending on the epoch interval used, 
the Spearman correlation coefficient is about 0.6-0.7. It is a 
stable correlation in the sense that the coefficient remains 
large when the epoch interval is perturbed. A correlation 
between light curve width and featu re 2 h as been previously 
discus sed by iBronder et al.l ()2008f ) and lArseniievic et al.l 
(12008ft . 

The pseudo-equivalent width for feature 2 also seems 
to be correlated with the fitted SALT colour parameter, 
as shown in the left panel of Figure 1121 but mainly in the 
epochs just after peak. While SNe with SALT-c > 0.2 do 
not appear to correlate with pEW, those below this rough 
limit seem to do. This could be interpreted as a sign of 
two different sources of reddening, where e.g. the highly 
reddened supernovae are dust extinctcd while most super- 
novae get their colour from some intrinsic property which 
is correlated with the strength of feature 2. 

As was shown in Sec tion 12. 1[ dust absorption according 
to ICardelli et a l. (1989j) would create a small pEW change 
in the opposite direction. 

To probe the origin of these lightcurve correlations the 
same epoch ranges and indicators as displayed in Figure [12] 
were examined using MLCS fit parameters. These results 
can be seen in Figure Q21 Correlation with lightcurve shape 
(A) is strong, while correlations with Ay are less clear. The 
correlation with A in the epoch range 0-8 is significant, 
while for the same epoch range, we find a correlation with 
SALT colour but only a weak correlation with stretch. 

The origin of these correlations is further discussed in 
Section [6j where we focus on feature 2, Si n A4000. 



shape (stretch) and spectral indicators seem to be present 
in the data analysed, we expect the host-lightcurve corre- 
lation to propagate to a correlation between spectral and 
host galaxy properties. Recent studies have also found in- 
dications of correlations between host galaxy properties, 
like mass and metallicity, and supernova absolute magni- 
tude that does n ot appear to be captured by lightcurve 
shape or colour ([Gallagher et al.l [20081: iKellv et all 120101: 



ISullivan et all 120101: iLampeitl et alj l2010bf ). It is thus of 
great interest to investigate if spectral indicators correlate 
with host galaxy properties, especially beyond what is re- 
lated to lightcurve stretch. 

All SN host galaxies were studied using the stellar 
formation code PEGASE. A descrip tion of this process 
can be found in ISmith et all (|201dh . and a comparison 
wi th lightcurve p r opertie s and Hubble diagram residuals 
in lLampeitl et al.1 (|2010b[ ). Here we use the estimated host 
galaxy type, host mass (in units of Mq) and specific star 
formation rate (sSFR; defined as the star formation rate 
per stellar mass, yr) and compare with spectral indica- 
tors. Host type is defined based on specific star formation 
rate: type zero indicate no star formation, type one mod- 
erately star forming ( —12.5 < \og(sSFR) < —9.5) and 
type two star forming (\og(sSFR) > —9.5). This param- 
eter thus largely overlaps with sSFR, with the important 
exception that hosts with no star formation (Type 0) are 
not displayed in the plots over sSFR values. 

As previously, this search for correlations was per- 
formed using various epoch intervals, but with only SDSS 
NTT/NOT SNe. Once again, the largest \r\ were found 
when probing the second feature (Si II A4000). When max- 
imising the correlation for host mass and specific SFR, we 
obtained the epoch interval between to 8 days past peak 
(which coincides with the epochs where a strong colour cor- 
relations is also seen). 

We also minimize the KS probability for indicator mea- 
surements from SNe in different host galaxy types to orig- 
inate from the same distribution. We do this through first 
compairing SNe from higly star-forming hosts (type 2) with 
remaining SNe (from host types and 1), and then com- 
pairing SNe from non star-forming hosts (type 0) with star- 
forming (type 1 and 2). 

The epoch interval with least probability of distribu- 
tions of pEWs from all host types being the same is epoch 
—9 to —2. These correlations are shown in Figure [HI 



5.2.1. Summary: Correlations with lightcurve parameters 

Pseudo-equivalent widths, as measured in this sample, do 
correlate with lightcurve properties. We recreate strong lin- 
ear correlations between f2-Si II A4000 both with SALT 
stretch and MLCS A. Correlations with SALT colour were 
also found. In general we see weak correlations between 
most pEWs around lightcurve max and both stretch and 
colour in the sense that wider, bluer SNe have small equiv- 
alent widths. 



5.3. Host-galaxy properties 

It is well documented that star forming, late type 
galaxies host brigh ter SNe with wider lightcurve shape 
(|Hamuv et aT1ll996b . Since correlations between lightcurve 



Before lightcurve peak The largest |r| connection between 
lightcurve width and pEW was found to be for f2 during 
the epochs right before lightcurve peak. In the left pan- 
els of Figure [14] we compare this feature with host galaxy 
properties. As can be expected, assuming a relation be- 
tween stretch and host galaxy type, we see that actively 
star forming galaxies have lower pEW f2 values than passive 
galaxies. A possible alternative explanation is that SNe in 
passive galaxies (Type 0) form a separate sub group: These 
all have large pEW f2 values, lower than average lightcurve 
widths and are only found in the very most massive host 
galaxies (as can be seen in the mid left panel of Figure [T4"|) . 
More statistics is needed to determine whether such a sub- 
group exists, or if a continuous trend with host type or mass 
is present. 
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Fig. 12. These correlation studies show ApEW, the measured pEW minus the expected pEW for the corresponding 
epoch (for the full sample of normal SNe la), vs SALT lightcurve parameter (stretch, colour) for different epoch intervals. 
These are given above each panel. The Spearman correlation coefficient is given in each plot. MC simulations show that 
the chance that these correlations are random consequences of a search through epoch ranges is one or less in 1000. 



After lightcurve peak We also found pEW f2 after lightcurve 
maximum to be related to host galaxy properties (right 
panels of Figure IT4|) . Lightcurve width and pEW f2 are 
correlated but not as strongly as before maximum. Instead 
we see a tentative correlation with SALT colour. 

Most significant in this epoch range is what seems 
to be a linear correlation between sSFR and pEW f2. 
Alternatively, as for the epochs before peak, this could be 
explained using a subgroup of SNe, in this case consisting 
of blue SNe with small pEW f2 values and high specific 
star formation rate. This relationship is emphasised if we 
include host mass information; all SNe in this group have 
low host masses. Note that all SNe in actively star forming 
hosts (Type 2) have smaller than average SALT-c lightcurve 
colours (c < 0.05), thus suggesting little dust extinction. 
Since we see a correlation with host galaxy mass, a random 
(uncorrelatcd) star formation rate would mean a correlated 
specific star formation rate (since this is the ratio between 
SFR and host mass). This clearly needs further study; we 
would expect SNe in small, star forming hosts to be more 
extincted. 

5.3.1. Summary: Correlations with host galaxy properties 

Host galaxy properties and spectral indicators, mainly 
pEW f2, are clearly connected. Of special interest is 
whether further subgroups among normal SNe can be iden- 
tified. Our results could be interpreted as a hint of the 
presence of two subgroups. One consisting of low stretch 
SNe with wide pEW f2, hosted by passive massive galaxies, 
and another consisting of blue SNe in actively star forming, 
low-mass galaxies. 

6. Discussion 

The discussion is split into a further examination of corre- 
lations in the Si II A4000 region fSection lG.ip . a search for 
signs of evolution with redshift (|6.2[) . a discussion of host 
galaxy properties f|6 . 3[) and finally we revisit some system- 
atic effects (|6.4p . This division does not mean that these 
topics are separate, they are rather closely related. 



We use composite (average) spectra as a tool for search- 
ing for physical differences between subsets of SN la spec- 
tra. Composite spectra can be normalised/created in a 
number of ways. This can make it hard to interpret the dif- 
ference between two composite spectra in terms of physics 
or compare composite spectra created using different mcth- 
ods0 

The composites shown here were created through scal- 
ing all spectra to have average fx = 1 in the studied region 
before combining. The uncertainties of the composites were 
then investigated through jackknife techniques. This sim- 
ple procedure is sufficient since our purpose is to compare 
different subsets of spectra, not provide a "true" supernova 
composite. 

6.1. Rest frame 4000-4500 

This region roughly correspond to features two and three. 
The absorption feature at ~ 4050, here called f2 and usu- 
ally attributed to Si II, has been shown to correlate with 
both luminosity and the stretch parameter in the sense of 
wider more luminous SNe having smaller equi valent width 
dBronder et al.l 120081: lArseniievic et all 12008ft . The same 
trend, using partially the same data, is seen clearly in this 
study using both SALT stretch (left top panel of Figure |T2j 
and MLCS A (left top panel of Figure We can further 
show that (i) the trend is continous with SALT-s/MLCS-A 
and (ii) strongest before lightcurve peak. 

Including the broader absorption region around 4200 A 
(here called f3 , attribu ted to Mgll amongst other ions), 
iGaravini etaH (|2007aD found a correlation between the 
"breaking point" of t his feature and lightcurve width and 
ISullivan et al.l (|2009|) found tentative signs of evolution, 
where high r e dshift objects have smaller equivalent widths. 
iBailev et all (|2009h used SNfactory data to look for the 
spectral regions most correlated with peak brightness; the 
best such was the ratio F(6420A)/F(4430A). 

5 For example, IVanden Berk et al.l (|200lh show how contin- 
uum measurements on Quasar composite spectra change de- 
pending on combination method. 
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Fig. 13. These correlation studies show ApEW, the measured pEW minus the expected pEW (for the corresponding 
epoch for the full sample of normal SNc la), vs MLCS2k2 lightcurvc parameter (A, Ay) for different epoch intervals. 
These are given above each panel. The Spearman correlation coefficient is given in each plot. The epoch ranges were 
chosen to correspond to those used in Figure [12j For the post lightcurve peak epoch range (0-8) we include correlations 
with both A and Ay- 



Finally one of the major spectral differences between 
normal and SN 1991bg-type SNe is the dominant Ti absorp- 
tion at these wavelengths, even though these subluminous 
SNe are not included in this analysis. 

This region is thus important both in order to under- 
stand and model supernova explosions as well as to make 
SNe better standard candles for use in cosmology. We will 
here try to probe the origin of the pre-peak correlation with 
stretch and the post-peak correlation with colour through 
comparisons of composite spectra. 

Pre-peak stretch composite All spectra contained in the 
epoch region between -6 and 1 (the region where the corre- 
lation with stretch was strongest) were included. Multiple 
spectra of single SNe were pre-combined. The sample was 
first divided into three stretch bins (low, intermediate and 
high). These composites can be seen in Figure [Tol Both the 
width and depth of the Si n A4000 feature grow smaller 
with increasing stretch (exemplified by the almost nonexis- 
tent absorption for overluminous SN 1991T-like objects). 
Around the epoch of peak luminosity, the photosphere 
temperature and composition seem to be such that the 
56 Ni yield, through the deposited energy, strongly affects 




Fig. 15. Composite spectra constructed from all peak spec- 
tra (NTT/NOT as well as reference SNe), divided according 
to stretch. The average lightcurve stretch, s, is printed in 
each panel, increasing from left to right. The average epoch, 
e, is also included for each composite. The plots show both 
the composite (thick red) and individual (grey) spectra. 
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Fig. 14. Pseudo equivalent width of feature 2 compared with host galaxy properties. Left panels: Pre peak epochs (-9 to 
2) with colours marking low/high lightcurve stretch SNe. Right panels: Post peak epochs (0 to 8) with colours marking 
low/high lightcurve colour SNe. Top: Host type vs pEW (0 not star forming, 1 moderately star forming, 2 highly star 
forming. Middle: Log of Host galaxy mass (M Q ) vs. pEW Bottom: Log specific start formation rate (log sSFR, j/r" 1 ) vs 
pEW. 



Si II A4000 ionisation and abundance, thereby causing a 
correlation between lightcurve width and Si II A4000 depth. 

As is dis cussed in | Sullivan et afl (|2009l ). changes in SN 
populations (jHowell et al.ll2007[) will be seen as a ch ange in 
average spectra with redshift. Bronder et al.l (|2008l) argues 
that this pEW-lightcurve width correlation could be used 



to correct SN la lightcurve instead of using the lightcurve 
width. Our results, using a larger sample, support this con- 
clusion. However, since photometry will be obtained for 
most SNe anyway and high S/N spectra in a small epoch 
range is needed, the practical use is limited. 
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Post-peak colour composite Correlating spectroscopic indi- 
cators with lightcurve SALT colour yielded a statistically 
significant correlation between colour and Si n A4000 pEW 
(using spectra obtained a few days after lightcurve peak). 
This correlation grows stronger if events with a colour above 
~ 0.2 are excluded. The need to remove highly reddened 
events could be explained if these are caused by a separate 
effect compared with moderately reddened SNe (e.g. cir- 
cumstellar absorption). A correlation between spectral in- 
dicators and colour would be of great interest, as this would 
show that at least some part of the colour-luminosity rela- 
tion seen in Type la SNe originates in the SNe and not in 
any extinction external to the explosion. 

In Figure [16] we explore composite spectra constructed 
based on colour. Excluding very reddened eventes (c > 0.3) 
the average SALT colour for spectra with rest frame epochs 
in the (0-8) range is 0.03. As with stretch one composite 
spectra was created out of all spectra with colour below 
average and one out of all spectra with colour above av- 
erage. The mean SALT colour values for the low colour 
composite is —0.03 and for the high colour composite 0.10. 

The most interesting difference in Figurc[T()]is a resolved 
small absorption at 4150A among the bluer SNe. This fea- 
ture can not be seen in the corresponding redder colour 
composite. Various authors have suggested absorption by 
C II, Cr II, Co II or Fe III in this region, e.g. dGaravini et al.l 
120051: iBranch et aLll2008t IScalzo et alJl2010h (as well as D. 
Sauer, private communication). The average pEW will be 
smaller for the blue composite since the peak at 4070 A pro- 
vides a bluer feature bound; this bound will be set ^4120 
A for the red (high colour) version. The average epochs for 
the two composites are very similar, and the effect is thus 
likely not created by epoch variations^ We finally note that 
the blue edge of the 3900 A peak is slightly weaker in the 
lower than average colour composite (left panel Figure [TBI) . 

Extinction by dust could, especially in the presence of 
noise, cause a pEW correlation with SALT-c or MLCS-Ay 
as follows: Highly extincted, i.e., highly reddened SNe, have 
a large fraction of the flux around 4000 A absorbed. The rel- 
ative S/N levels will thus be smaller for redder SNe, which 
could make it impossible to resolve the small 4150 A ab- 
sorption. This effect would be more significant for SNe with 
already low S/N. 

While it is clear that the composites differ with respect 
to the 4150 A absorption, we cannot further distinguish 
between: (i) random small sample effects (ii) lower S/N 
preventing resolution of the small dip and (iii) a difference 
in the nature of the explosion correlated with SALT colour. 
The last alternative is intriguing, and strengthens the in- 
terest in using spectral indicators to disentangle the depen- 
dence of supernova properties on colour. 

6.2. Evolution 

Evolution of SNe la is a "thorny" issue; it is hard to 
make firm predictions or document that SN properties 
have "evolved". From a theoretical point of view, chang- 
ing metallicity (and other parameters changing with time, 
like galaxy ages) could have an impact on SN la pro- 
genitor systems, but the details remain largely unknown. 




pEW f3 = 16 A <e> = 4.83 

<c> = -0.034 



4000 4200 
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Fig. 16. Composite spectra constructed from all spectra 
with epochs in the range (3-8). The left panel composite 
is based on spectra with colour below average, the right 
panel one is based on spectra with colour above average. 
The plots show both composite (red) and individual (grey) 
spectra. The average SALT colours and epochs are printed 
for each composite, as well as the pEW for the indicated 
region. The dashed square marks a region where a small 
feature seems to exist only in the low colur sample (SALT 
c < 0.03), either because it is not resolved or because it 
is very rare among the higher colour SNe (0.03 < SALT c 
< 0.3). 



Previous work on SN la indicates that their properties ap- 
pear to vary depending on the surrounding environment 



Feature 3 also show some correlation with colour in the same 
epoch range. Since these features share a border this correlation 
can very well be caused by flux change around 4100 A. 



pear to vary dependi ng on tne sur roundin g environment 
(iGallagher et al.ll20q8t ISullivan et alJIgOldl: iLampeitl et al.l 
12010b! ) and redshift ()Howell et alj2007l) . So far these effects 
are interpreted as due to population evolution (a larger frac- 
tion of "prompt" hig h stretch objects at higher redshifts, 
iMannucci et alj|2005h . 

Such evolution would not per se imply any fundamental 
problems for SN cosmology. However, if spectral evolution 
not captured in lightcurve shape fits would be detected, 
this could create a cosmological parameter bias since these 
spectral changes are likely to affect luminosity. We thus 
need to show that we understand such changes. 

Both the ESSENCE dFolev et al.l l2008h and SNLS 
(|Balland et al.l 120091: ISullivan et al.l l2009t ) collaborations 
have reported indications of evolution, or systematic dif- 
ferences in redshift binned composite spectra. Both groups 
use samples where the average lightcurve width increases 
with redshift. Based on the lightcurve shape correlations 
found above, we would thus expect decreasing equivalent 
widths with redshift. This agrees qualitatively with the dif- 
ferences found by both groups. 

6.2.1. The pEW-deficit subsample 

As was discussed in Section 14.21 the deviations between 
the reference and NTT/NOT SNe can be attributed to a 
set of SNe with lower pseudo-equivalent widths than the 
non-peculiar reference SNe. These were collected in a pEW- 
deficit sample. We will here study why these NTT/NOT 
SNe deviate and whether this is a sign for SN la evolution. 
As our basic tool we create composite spectra of deficit 
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SNe and compare with composite spectra based on normal 
non-deficit NTT/NOT SNe. 

As can be seen in Figure [51 most spectra with low pEW 
f2 can be found before lightcurve peak while the deviations 
in f4 are most clear after peak. We therefore make two sep- 
arate studies: In Figure [171 we compare spectra with deficit 
pEW f2 in the epoch range —5 to with non-f2deficit spec- 
tra in this epoch range. In Figure [18] we compare spectra 
with deficit pEW f4 in the epoch range to 5 with non- 
f ^deficit spectraQ In each figure we also include the mean 
and dispersion of a number of properties for each subset. 
These include both lightcurve properties and estimates for 
possibly systematic effects like slit loss and host galaxy con- 
tamination (See lOstman et al.l (|2010l ) regarding these esti- 
mates). We caution that most composites consist of compa- 
rably few objects (<~ 10), and are thus sensitive to random 
fluctuations. 
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Fig. 17. Pre peak, f2 deficit vs normal: Composite spec- 
tra of "reference normal" SNe (dashed blue) compared with 
pEW-deficit SNe (solid red) in the epoch range —5 to days 
past peak brightness. The shaded regions correspond to 
jackknifc dispersion (light blue for normal, grey for deficit). 
The do tted orange line corre sponds to the deficit compos- 
ite with lCardelli et all (jl989t) dust applied that matches the 
colour difference. The arrows indicate features discussed in 
the text. Mean properties of the composites are stated in 
the Figure. 



Description of (physical) differences The mean spectra are 
similar in both cases, except for some limited regions. These 
are marked with black arrows. Starting with Figure [T7] (the 
pEW f2 deficit before peak brightness) the deficit com- 
posite has more flux (less absorption) at the silicon lines 
Si II A4000 (f2) and Si n A5800 (f5), as well as at 4400 A. 

In Figure HH (pEW f4 deficit after peak brightness) we 
find more flux (less absorption) in the normal composite 
both at 4150 A and 4500 A. 



7 Spectra of three SNe (17880, 16215 and 18466) were too 
deformed (reddening/slit loss) to be used in a composite and 
were removed from the analysis. 
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Fig. 18. Post peak, deficit vs normal: Composite spec- 
tra of "reference normal" SNe (dashed blue) compared with 
pEW-deficit SNe (solid red) in the epoch range to 5 days 
past peak brightness. The shaded regions correspond to 
jackknife dispersion (light blue for normal, grey for deficit). 
The dotted orange l ine co rresponds to the deficit composite 
with lCardelli et ail ((1989f ) dust applied matching the colour 
difference. The arrows indicate features commented on in 
the text. Mean properties of the composites are stated in 
Figure. 



Do the differences agree with lightcurve correlations? We can 

understand most of the observed differences seen based 
on the correlations between pseudo-Equivalent widths and 
lightcurve properties. The pEW f2 deficit SNe in Figure [T71 
have, on average, higher stretch (1.03) than the normal 
set (0.92). As was seen above, pEW f2/Si n A4000 has a 
strong negative cor relation with stretch . The same is true 
for Si II A5800 fsee iNueent et al.lll995l) . Interestingly, the 
region around 4400 A show a similar variation as the Si 
lines. This r e gion w as one of the flux ratio "legs" used by 
iBailev et al.l ( 20091 ) to standardize SNe la (as an alternative 
to lightcurve width). This suggests that the same physical 
process is responsible for the absorption seen here as at the 
Si lines. 

The f4 composites in Figure [TS] have a (small) difference 
in average stretch, but also in SALT colour: The average 
colour for the deficit sample is —0.02 while it is 0.03 for the 
normal sample. The differences in colour is smaller than the 
dispersion, but does agree in direction with what we expect 
form the tentative colour correlation seen in Tabled As the 
NTT/NOT SNe have a lower mean colour this could explain 
why SNe at higher redshifts seem to have smaller pEW f40 

We also find a significant difference at 4130 A, possibly 
related to the discussion (also including non-SDSS spectra) 
in Section IQl and Figure [TBI More flux is absorbed among 
low-reddened objects. 



8 We can make a ver y speculative extended comparison with 
the IBailev et all (|2009t ) results: If we assume that lightcurve 
width causes spectral differences as seen in Figure [17] and red- 
dening causes spectral differences as seen in Figure \TE\ and we 
want to find one wavelength where flux correlate both with 
lightcurve width and r eddening, we are le d to the small over- 
lap region at ~4400 A. IBailev et all (|2009f ) see most significant 
correlation at 4430 A. 
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We finally note that some pEW selected subsamples 
have a smaller absolute magnitude dispersion (restframe 
B band). Sample sizes are too small to properly evaluate 
this effect. 



Could the differences be caused by systematic effects? 

Systematic effects such as inaccurate correction of host 
galaxy contamination, slit loss or dust extinction could sig- 
nificantly affect spectral comparisons. 

To study dust extinction we have applied [Cardelli et al.l 



(|1989l) type dust matching the difference in SALT colour to 
the (less reddened) deficit composite (dotted orange line in 
Figures [T71 and ITS)) . In none of the cases does this approach 
create a better match to the normal set. Standard dust 
could thus not by itself create the observed differences. 

In each figure we have also included mean and disper- 
sion for our (conservatively) estimated slit loss and host 
contamination for each subset. Slit loss is defined as the 
fractional flux loss at 4000 A (observed frame), host galaxy 
contamination as the fraction of galaxy flux in the observed 
g band before host galaxy subtraction. Slit loss values are 
consistent among all subsets. Host contamination levels arc 
almost identical for the pEW f2 samples shown in Figure H71 

For the pEW f4 divided subsets in Figure [TS] the mean 
contamination level is lower for the deficit sample (~ 10%), 
but it is still comparably low for the reference SNc (~ 
19)%). It is very unlikely that this would cause a host con- 
tamination systematic difference since (i) the spectral com- 
posite differences seen are localized to small wavelength re- 
gions and (ii) it is the less contaminated objects that seem 
to deviate. 

In summary, we find no reasonable combination of sys- 
tematic effects that could create the deficit subsample. 



Are the pEW deficit SNe "peculiar"? The deficit NTT/NOT 
SNe were chosen since they do not correspond to normal 
reference SNe la. It is natural to ask whether they instead 
correspond to SNe la locally defined as peculiar. None of 
them seem to be "very" odd, in the sense of SN 2000cx 
or SN 1991T. But they do correspond to several other 
slightly peculiar Shallow Silicon (SS) SNe, like SN 19 99aw, 
SN 1999bp and SN 1999bn (See iBranch etall (120081) for a 
definition of SS). The deficit spectra have, in general, wide 
lightcurves and small "shallow" silicon features. 

In total we have 21 SNe belonging to either the f2 or the 
f4 deficit subgroups. Among these we find three NTT/NOT 
SNe that can be classified as SS SNe (15132, 17497 and 
19899). Seven (out of 21) SNc have spectra good enough 
to be classified as not of this subtype, leaving 14 SNe as 
possibly peculiar. We thus have between 3 and 14 SNe that 
would have been identified as peculiar if observed locally. In 
total (deficit and normal) we have 41 spectra in this epoch 
range, which translates into a fraction of SS SNe b e tween 
7% (3 out of 41) and 32% (14 out of 41). iLi et al.1 (pOlOh 
found that 18% of all SNe in a luminosity limited sample 
was of the SN1991T subtype but that it gets significantly 
harder to identify these without early spectra. 



Summary: Evolution The "deviations" between the 
NTT/NOT data set and the normal reference set can 
be explained through a combination of two (connected) 
effects: 



A fraction of "borderline" peculiar SNe, mainly similar 
to Shallow Silicon SNe, that would have been identified 
as such if obscreved locally and thus do not exist in the 
reference sample of normal SNe. 

SN features change with lightcurve parameters. Sets 
with different lightcurve parameters will thus have dif- 
ferent spectral properties. 



6.2.2. Sensitivity to evolution models 

While we cannot exactly determine our sensitivity to (an 
unknown) spectral evolution, we have simulated how well 
we would have detected some models of evolution. This pro- 
cess is further described in Appendix IA"1 These simulations 
used all property distributions in the NTT/NOT sample 
and determined how well we would detect changing SNe 
subtype distributions, that is if the fraction of evolved SN 
increased with redshift. We conclude that most of the mod- 
els studied should have been discovered at least at low sig- 
nificance (<~ 2a), using one or more indicators and possibly 
removing high bias events. These models are, however, not 
realistic. It is possible that an increased fraction of "deficit" 
SNe is a sign of evolving subtype distribution, but if so this 
does not limit the use of SNe la as standard candles (since 
these SNe seem to follow the same luminosity-width rela- 
tion as other SNe la). 



6.3. Host galaxy properties 

Different galaxy types give rise to different Type la SN 
populations. First, the distribution o f lightcurve param- 
eters differ. Se cond, as i ndicated by Kelly et al.l (|2010t ): 
ISullivan et ail (|2010f ) and iLampeitl et al.l (|2010bD . Hubble 
diagram residuals seem to correlate with host types. 
However, it is still not clear how strong these effects are 
or the causes. 

The comparisons between Si II A4000 and host galaxy 
properties could help clarify these questions. As is discussed 
in Section [53] we see indications of relations between spec- 
tral indicators and host galaxy properties. 

For the epochs after lightcurve peak, we see signs of 
correlation with both host mass and specific star forma- 
tion rate. The measurements could also be interpreted as 
belonging to different subtypes of SNc la (e.g. a subset of 
SNe produced in low mass hosts with high sSFR whose 
spectra after peak show very small Si II A4000 values). 

Correlations with host mass could be present also for 
pre-max spectral epochs. This is best seen by separating 
low and high lightcurve stretch SNe. 

We have also compared the host galaxy types with the 
subset of SNe defined as (possible) shallow silicon (SS) SNe. 
We find that all clear SS SNe originate in actively star form- 
ing galaxies (Type 2) as well as three out of four likely SS 
SNe. These six SNe constitute almost half of all the 13 SNe 
with host galaxies of Type 2 (with spectra in the —5 to 5 
epoch range). Considering the SNc too noisy to be iden- 
tified as belonging to a subtype this implies that half or 
more of all SNe in our sample from actively star forming 
galaxies are similar to Shallow Silicon SNe. We also note 
that these SNe tend to have blue SALT lightcurve colours 
and originate in lower mass host galaxies. 
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6.4. Systematic effects 

We have studied several possible systematic effects. In 
Figure [TO] we show some sample indicator measurements 
with NTT/NOT SNe divided into subsamples according to 
possibly systematic effects. For example, if host galaxy con- 
tamination would yield a systematic bias we would expect 
low contaminated SNe differ from the high contamination 
set. We did not detect any such differences or signs of sys- 
tematic differences. 

Use of local templates Wc make use of our knowledge of 
local SNe and SN templates in a number of different ways. 
This includes both the host-galaxy subtraction method 
used here and identification of SNe using templates (e.g. 
using software like SNID). 

This approach is clearly not ideal and creates tension 
when probing for evolution with redshift. A few items to 
note regarding this: 

— For moderate evolution, local templates should provide 
a fair match, and we should be able to accurately type 
SNe and extract clean spectra. It would be a striking co- 
incidence if the evolution could be completely masked 
through varying the three galaxy eigenspectra used to 
model the galaxy. The coefficients of the fit are strongly 
ovcr-dctcrmincd by the number of wavelength bins fit- 
ted. 

— It is significantly harder to identify completely new sub- 
types among low S/N distant SNe. While this would be 
very interesting it is not within the scope of this study. 

— Most cosmological surveys use an identification mecha- 
nism relying on known templates to select the SNe to 
include in samples for cosmological fits. 

— The synthetic spectra discussed in Appendix [A] are con- 
structed from SN spectra which differ from the tem- 
plates used in the host subtraction. Simulations show 
that we can still get correct indicator measurements af- 
ter subtractions. 

We thus conclude that the use of SN templates is not a 
fundamental objection when searching for evolution among 
SNe la used for cosmological studies. Our general conclu- 
sions would not have changed if only objects with small 
contamination were included, see Figure 1191 although they 
would have had smaller statistical significance. 

Noise and slit loss Two further possible systematic effects 
are noise and slit loss. The possible bias effects when com- 
paring high and low S/N data should not be underesti- 
mated. To probe these effects, the NTT/NOT sample was 
split into high/low S/N and high/low slit loss samples and 
the results compared. In Figure Q1J] wc show these for fea- 
ture 3. No significant bias was detected. Simulation results 
show that we can adjust indicator errors according to noise 
levels. 

However, we note that noise can still limit our abil- 
ity to determine whether spectral differences are "real" . 
In Section 16.11 we discussed whether higher noise levels in 
more reddened SNe could hide small spectral features seen 
among less reddened SNe. 

Selection effects The NTT/NOT sample was obtained as 
part of the SDSS-II Supernova Survey, a rolling search. Very 



faint targets were usually scheduled for typing at larger 
aperture telescopes. Even though the NTT/NOT sample 
has a normal distribution of lightcurve parameters, it is thus 
possible that these SNe arc not representablc for the SN la 
population as a whole (e.g. missing some faint objects). We 
have not attempted a full study of the completeness of the 
NTT/NOT sample. 

The reference sample is very inhomogencous. These SNe 
were not observed as part of a rolling survey, and more 
or less peculiar objects are over represented. Also, since 
multiple spectra exist of many objects, these objects will 
carry larger weights. 

Finally, the distribution of epoch and colour differ be- 
tween the samples. We can thus not expect full agreement 
between the samples. 

As we conclude above the differences between the refer- 
ence and NTT/NOT sample can be explained by a fraction 
of semi-peculiar Shallow Silicon SNe. To properly identify 
these objects, multiple high S/N spectra at early epochs are 
needed; it is thus not surprising that these are identified in 
low-z data but not in the NTT/NOT sample. 

7. Conclusions 

We have measured both pseudo-Equivalent widths and line 
velocities of individual optical spectra observed at the NTT 
and NOT as part of the SDSS-II Supernova Survey. These 
spectra cover the redshift range 0.05 — 0.3. Our spectra were 
compared with a low-redshift sample to probe a possible 
evolution between local SNe and SNe at cosmological red- 
shifts. The samples were then combined and all SNe were 
used to investigate possible correlations with lightcurve 
properties. 

The differences between reference and moderate redshift 
SNe can be well described by a fraction (<~ 20%) of slightly 
peculiar, possibly Shallow Silicon, SNe la. 

The linear correlation between Si n A4000 pseudo 
Equivalent-width and lightcurve shape is very significant, 
both when using SALT stretch and MLCS A parameter- 
isation. We also found correlations between this feature 
and SALT lightcurve colour (particularly if highly reddened 
events are excluded) in spectra observed roughly during the 
first week after lightcurve peak. This could be an effect from 
intrinsic colour dependence or a sign of different noise lev- 
els. In this epoch range, Si n A4000 correlates with MLCS 
A, but we also see a faint correlation with MLCS Ay for 
NTT/NOT SNe. Further studies have to conclude whether 
an intrinsic physical correlation with reddening exist. 

We do not see any significant correlation between spec- 
tral properties and absolute magnitude, but we do find a 
smaller magnitude dispersion among SNe subsamples de- 
fined through pEW. 

We also found connections between host galaxy prop- 
erties and spectral indicators. As for the correlations with 
lightcurve parameters, these seem strongest for Si II A4000. 
Future studies are needed to confirm whether these are real, 
and whether they, in turn, derive from lightcurve parame- 
ters. The correlations could be explained by a subset of SNe 
with weak Si II A4000 appearing in hosts with low mass and 
high star formation. 
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Fig. 19. Study of systematic effects: Each panel highlight pEW measurements possibly affected by systematic effects. 
Top left: Feature 3 highlighting low contamination events. Top right: Feature 4 highlighting low contamination events. 
Bottom left: Feature 3 highlighting high S/N SNe. Bottom right: Feature 3 highlighting low slit loss SNe. No major 
deviations in subsets with possible systematic errors are detected. 
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Appendix A: Host-galaxy subtraction uncertainties 
and evolution detection limits 

Estimation of the SED of contaminating host galaxy light is 
an essential step if spectral indicators in contaminated and 
uncontaminated spectra are to be compared. This will, in 
turn, be unavoidable when comparing nearby (usually with 
the SN clearly separated from the host galaxy core) and 
distant SNe (where SN and galaxy light are degenerate). 
In lOstman et al.1 (|2010t ) we present the host-galaxy sub- 
traction pipeline applied to the NTT/NOT SNe. In short 
this method consists in matching a SN template with a 
number of galaxy eigencomponent spectra, including a slit 
loss/reddening correction. Even if the observed SN SED de- 
viates slightly from the SN templates used in the fit, the 
very large number of wavelength bins compared with the 
few fit parameters (five) will allow this SN deviations to 
remain after the subtraction. 

However, the host subtraction produces an increased 
indicator measurement uncertainty and possibly a bias. It 
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is important that this uncertainty or bias is estimated. In 
this Appendix we describe the extensive simulations that 
were run to study the effectiveness of the host subtraction. 
These simulations were used to calculate a systematic bias 
and uncertainty for every measurement, depending on the 
shape of the indicator and contamination level. 

As a second step of these simulations we used suggested 
(metallicity) evolution models to study under which cir- 
cumstances these would be detected assuming the proper- 
ties of the NTT/NOT data set. 

A.l. The subtraction pipeline 

The subtraction pipeline is described in detail in 
lOstman et al.l (|2010l ). The input parameters are flux den- 
sity and (optionally) error, observer frame wavelength, red- 
shift and an epoch estimate. This pipeline thus operates 
identically for real and simulated spectra. A range of in- 
ternal fit parameters can be changed, including which tem- 
plates and host galaxy eigencomponcnt spectra are used 
as well as the nature of slit loss/extinction approximation. 
The fit parameters were optimised and fixed during a series 
of test runs. 

A. 2. Synthetic spectrum simulations 

To test the reliability of the estimated host galaxy spectra 
and the impact on spectral indicators, a large number of 
simulated contaminated spectra were created. Besides con- 
tamination, these simulations included realistic slit loss and 
noise levels. The synthetic spectra are created from 

— a supernova spectrum 

The SN spectra used as templates all have high S/N and 
low contamination. Their epochs arc similar to the ones 
of the NTT/NOT spectraQ Eleven different SN spectra 
are used: five of SN 20 03du (epochs -6, -2 4, 9, 10, 17) 
(IStariishev et al.ll2007t ). one of SN 1998aq (jBranch et al.l 
l2003h at peak brightness, two of the subluminou s SN 
1999by (epochs -5 and 3) (jGarnavich et al.ll2004l) and 
two of t he peculiar and lumi nous SN 1999aa (epochs -5 
and 0) (iGaravini et al.ll2004Tl . 
reddening is added to the SN spectrum 
The reddening is added using the lCardelli elaTl (fl989h 
extinction law using a total-to-selective extinction ratio 
Ry of 2.1 and a colour excess E(B — V) drawn from the 
distribution of E(B - V) obtained from the NTT/NOT 
lightcurve fits, 
a galaxy spectrum 

Four galaxy te mplates of varying ty pe (elliptical, SO, 
Sa and Sb) from lKinnev et alT(|l996t) are used together 
with three real galaxy spectra observed at NTT at the 
same time as the SN spectra analysed here (host galaxy 
spectra for SDSS SN7527, SN13840 and SN15381). The 
contamination level is randomly chosen between and 
70% for the g band. These simulations were later ex- 
tended in a second series where 50 randomly chosen 
SDSS galaxy spectra were used. Figures displayed here 
are based on the first run series, but results are similar 
when including the second set of galaxy spectra. 



These templates are omitted from the subtraction pipeline 
since this would make the fit trivial. 



redshift 

The object redshift is randomly drawn from the 

NTT/NOT redshift distribution. 

slit loss is added to the SN spectrum 

The different ial slit loss functions are taken from 

lOstman et al.l (|2010t ) and correspond to typical 

NTT/NOT situations and range from insignificant to 

severe. 

noise addition 

A S/N value is randomly chosen from the NTT/NOT 
spectral S /N distribution. Poisson noise is added to the 
spectra until the target S/N is achieved. The shape of 
the noise is determined as a linear combination of the 
input spectrum and a randomly chosen NTT/NOT sky 
spectrum. The linear combination is regulated such that 
the highest S/N value in the NTT/NOT sample corre- 
sponds to no contribution from sky noise, the lowest 
S/N corresponds to complete dominance by sky noise 
and intermediate values to a combination of the two 
error sources. 



All of the created synthetic spectra were then processed 
through the host subtraction pipeline and spectral indica- 
tors were measured. The measured spectral indicators could 
then be compared with the ones obtained from the original 
SN spectrum. The subtractions were thus evaluated only 
with respect to how well correct indicators were measured. 



A. 3. Simulation results 

The simulation results can be analysed in a number of ways: 
Looking at specific SN spectra, specific galaxy types, spec- 
tra with more or less slit loss or contamination or any com- 
bination of these. For each of these subgroups errors in all 
equivalent widths and velocities can be calculated. 

In general simulations are stable with the following char- 
acteristics: A small bias for very low contamination lev- 
els that decrease with added contamination and a random 
dispersion that increases with contamination. The size of 
these effects vary slightly from feature to feature. The small 
bias for low level contamination means that the subtraction 
pipeline finds "something" to subtract even when no con- 
tamination was added. This is fully consistent with having 
a small amount of host light already present in the template 
spectra. But we cannot rule out that that a part of this bias 
is caused by the subtraction methods. In practise we do not 
perform host subtraction on spectra with very low contam- 
ination levels. In all cases the full bias as estimated in the 
simulations is retained, thus generally overestimating the 
bias levels. 

Sample simulation results are presented in Figure IA.1I 

The simulations were evaluated with and without added 
noise. Noise was found to increase the error dispersion but 
not introduce any significant bias. The added dispersion 
was comparable to uncertainties yielded from the desig- 
nated noise simulations. We thus separate errors from host 
contamination and noise. See Appendix iBl for a further dis- 
cussion about noise and filtering. For final spectra the sys- 
tematic uncertainties will be the sum in quadrature of the 
respective subtraction and noise systematic uncertainties. 
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Contamination 9% 
Avg error 5% 
Prms error 11% 



-0.2 0.2 

Fractional error 



Contamination 26% 
Avg error 3% 
Prms error 12% 



-0.2 0.2 

Fractional error 



Contamination 44% 
Avg error 1% 
Prms error 13% 



-0.2 

Fractional error 



Contamination 61% 
Avg error -2% 
Prms error 15% 



-0.2 0.2 

Fractional error 



Fig. A.l. Sample host contamination simulations results. For every simulated spectrum the final fractional error is 
calculated (fractional error is used so that all templates of different epoch and subtypes can be added and analysed 
as function of contamination). The panels show the distribution of errors, divided into four contamination bins (0 — 
17.5,17.5 — 35,35 — 52.5 and 52.5 — 70% in g band). The average contamination, average error and Population RMS 
(Prms) is printed for each bin. The average error (shown as dashed orange line) indicates a small bias, decreasing with 
contamination. The dispersion indicates a random error from host subtraction, increasing with contamination. These 
plots are based on pEW f3; other pEWs show similar results. 



A. 4. Alternative subtraction methods 

A number of alternative host subtraction methods were 
tried. These included two fundamentally different fitting 
methods: Linear fits using all nearby SN spectra as SN 
templates and photometry fixed galaxy subtraction where 
the host galaxy photometry is used to constrain the galaxy 
shape and proportion. Both methods relax the dependence 
on the SN template, the first through including a larger 
variety of such and the second through not using any tem- 
plate at all. However, in general the multiplicative method 
including the slit loss/reddening correction was found to be 
superior in most cases and generally more stable. 

A number of different implementations of the subtrac- 
tion pipeline were also tried. These included modifying the 
number of galaxy eigenspectra, the origin of these eigen- 
spectra and changing constraints on the eigenspectra pro- 
portions. The host galaxy subtraction method described 
above was the final product of these tests. 

However, there will be individual objects, for which the 
host subtraction fails or performs less than ideal. This is 
a natural consequence of the degeneracy between SN, host 
galaxy and noise. For some of these objects alternative sub- 
traction methods could have been better suited, but for 
consistency uniform host subtractions were used. The sim- 



ulations were designed to estimate the bias caused by such 
failed subtractions. 

A. 5. Evolution models 

Since it is unknown if evolution exists and how it, if existing, 
affects the SN la SED, it is impossible to predict whether 
evolution could be detected with the NTT/NOT SNc. But 
we can still study proposed models to quantify how well 
these effects would be detected. Two different models were 
considered here: First ad hoc decrease of the depth of fea- 
ture 3 and 4, where the frac parameter regulates the per- 
cent decrease of these depths. This modification was in- 
spired by the indic ation of changes in these fe atures seen by 
iFolev et all (|2008l) and ISullivan et all (|2009l) . As a second 
set of models we use the spectral changes cause d by one low 
and o ne high metallicity model simulated by iLentz et al.l 
(|2000D . For spectrum templates with epochs less than —2.5 
the 15 days after explosion model was used, otherwise the 
day +20 model. 

All base SN templates used in the above simulations 
were modified according to the evolution models, and pro- 
cessed through the subtraction and measurement pipelines 
again. The modifications as applied to the SN spectrum 
of SN2003du observed at April 30 2003 is displayed in 
Figure IA.2I 
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These models should not be considered realistic evolu- 
tionary models to be tested, but rather tests as to what level 
of evolution can be detected assuming host subtraction un- 
certainties. They are however, examples of evolution that 
would not be detected by visual inspection of noisy data 
but could still effect SN la cosmology. 

A. 6. Evolution detection limits 

All measurements on "evolved" host galaxy subtracted 
spectra are collected and compared to the true unevolved 
reference values. This difference between measurements can 
then be compared with the estimated statistical and sys- 
tematic uncertainties and the likelihood of detecting the 
evolution studied. Sample evolution detection probabilities 
for evolved SNe is shown in Figure IA.3I 

These comparisons show that most evolved SNe would 
be detected. However, the detection limits we are searching 
for must be realistic: We do not expect all SNe at higher 
redshift to be evolved, but rather the fraction of e.g. low 
metallicity SNe will change. To study this limit we designed 
a further simulation based on the NTT/NOT redshift dis- 
tribution. The probability of each SN to be evolved accord- 
ing to one of the above models, is set to be proportional 
to redshift and reach 50% at the average redshift of the 
NTT/NOT data set. For each model we repeat the mea- 
surement 5000 times and in each we randomly select which 
SNe are evolved. The total spectral indicator offset is cal- 
culated and compared to the uncertainties, thus obtaining 
a distribution of the evolution detection limit. 

In Table ED detection limits assuming all NTT/NOT 
SNe (including high contamination) are listed for a num- 
ber of indicators for the models for evolution/metallicity 
discussed above. These limits are completely dominated by 
the systematic bias levels of the high contamination events, 
since the systematic bias is set to be a systematic floor 
where the largest bias contained is used. A more realistic 
and less conservative estimate arises when we remove the 
highest bias/contamination events; these limits are given in 
Tabic [O 

These results show that we would be sensitive to all but 
the very weakest of these evolution models using at least 
one indicator, albeit at a fairly low significance level. 

Table A.l. Probability of detecting models for SN la evo- 
lution. Each column corresponds to one model (first col- 
umn is no evolution) , see text for further description. Each 
row corresponds to a search for evolution using the spec- 
ified spectral indicator, assuming the population changes 
linearly with redshift. Numbers are the detection level in 
standard deviations using max statistical/systematic error 
as irreducible global error . The last line is an example 
where measurements of two indicators are combined to in- 
crease sensitivity ('*' = comparison not made). 
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Fig. A. 2. Models of evolution/metallicity changes applied 
to SN2003du. f (frac) models consist of a decrease in the 
depth of th e f3 an d f4 features, Z 0.5 corresponds to the 
iLentz et al.l (|2000D model of increase d metallicity, Z -1.5 
corresponds to the ILentz et al.l (|2000l ) model of decreased 
metallicity. 

Table A. 2. Probability of detecting models for SN la evo- 
lution. Each column corresponds to one model (first col- 
umn is no evolution), see text for further description. Each 
row corresponds to a search for evolution using the spec- 
ified spectral indicator, assuming the population changes 
linearly with redshift. Numbers are the detection level in 
standard deviations when removing highest bias events. The 
last line is an example where measurements of two indica- 
tors are combined to increase sensitivity. 
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A. 7. Velocity host subtraction errors 

Host contamination could affect velocity measurements ei- 
ther through introducing a false minimum or through mod- 
ifying the position of the true minimum. Studies of simu- 
lated spectra show that velocity errors do increase with 
contamination, but below an r-band contamination of 60%, 
the errors are small compared to statistical and noise un- 
certainties. 

Host subtraction methods in general perform similarly. 
The same subtractions as for pEWs are used (for consis- 
tency) . Systematic uncertainties as estimated from the sim- 
ulations are added to all measurements. 

Appendix B: Filtering and uncertainties due to 
noise 

Random noise will degrade data quality, making measure- 
ments less accurate. For low S/N SN spectra, the con- 
ventional solution is to apply a filter to remove the high- 
frequency noise. This technique works well if small levels of 
filtering are used (filtering/smoothing are considered iden- 
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Indicator difference compared with no evolution (%) Indicator difference compared with no evolution (%) 

Fig. A. 3. Sample study of how well evolution is detected in simulated spectra. The "30 %" evolution model was applied 
to all template spectra and the measured indicators compared with the unevolved measurements. The panels show the 
distribution of fractional difference, divided into the same contamination bins as in Figure lA. 11 The total uncertainty in 
each bin (bias and dispersion) as estimated above is shown as an orange dashed line. Events where the reported difference 
is larger than uncertainties would be seen as deviating. In this sense detectable events are shown as hashed bins. The 
fraction of detected events is shown in each panel. This fraction decreases with contamination. 



tical processes here), where the true shape is clearly visible. 
For noisy data it is no longer obvious what filter to use or 
how accurate results are. 

According to the definition, pseudo-equivalent widths 
run from one wavelength extremum point to another. This 
makes such measurements extremely sensitive to noise: if 
any noise peaks remain, the pseudo continuum will be de- 
fined from there. To remove these, and create unbiased 
data, strong filtering is needed for low S /N data. We would, 
however, not want to filter high S/N spectra (at any red- 
shift) too much since this would destroy information. We 
would also like to estimate noise uncertainties. 

A further complication caused by filtering is that errors 
in filtered bins are correlated. 

A series of Monte Carlo simulations were run in order to 
(i) compare filter methods, (ii) determine filter parameters 
and (iii) estimate associated uncertainties (while avoiding 
having to determine filtered error correlations) . These sim- 
ulations are described below. 

B.l. Filter method comparison 

Three filters easy to implement are (1) the boxcar filter, 
which is simple averaging over a wavelength range, (2) the 
variance- weighted Gaussian filter where the smoothed value 



in a pixel is determined from a surrounding region weighted 
by a Gaussian determined by the inverse variancq3 and 
(3) the FFT filter, where all frequencies above a certain 
maximum frequency are removed from the spectrum. 

In order to determine which filter method works best 
and find optimal filter parameters, MC simulations were 
run. Random noise was added to template SN spectra af- 
ter which the S /N was determined, the spectra filtered and 
indicators measured. For each method the optimal filter 
parameters were found through minimisation vs. the true 
value. This process was repeated until MC errors were suffi- 
ciently small. It was found that there is no optimal method 
with a single set of parameters that worked over the com- 
plete range of varying features and S/N values. All methods 
can yield non-biased values if correct filter parameters are 
used. The correct filter parameters should be determined by 
the actual noise level and the nature of the feature studied 
(broad or sharp). 

Since all methods can be made to work but none will 
work with a single set of parameters, we selected the sim- 
plest method, the boxcar filter, as described below. 



See lBIondfn et all (|2006T ) for a more detailed description. 
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Filter MC: pEW error f3 epoch 



Filter MC: pEW error f7 epoch 
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boxcar width (A) 



20 40 60 
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Fig. B.l. Average pEW error for feature 3 (left panel) and feature 7 (right panel) at the epoch of maximum light. The 
noise level, expressed through the logarithm of the pseudo-S/N, increases along the y-axis and the filter strength along 
the x-axis. Darker shades show smaller errors. 



B.2. Optimal boxcar filter parameters for pEW measurements 

The above simulations showed that true pseudo-equivalent 
widths can be measured from noisy spectra after binning, 
but correct bin widths must be used. A range of MC sim- 
ulations were run to determine the widths to use and the 
typical errors caused by noise. This procedure is detailed 
below. 

Noise was generated with a certain amplitude. A grad- 
ually stronger filter was applied, while measuring relevant 
features at each stage. Through comparison with the true, 
noiseless values, the errors are obtained. For each iteration 
a "pseudo-S/N" is calculated as follows: A minimal boxcar 
(spanning three bins) is applied, and a pseudo-S/N can be 
calculated by comparing this with the original spectrum. 
This value serves as an initial estimation of noise level, and 
can later be compared to real spectra (adjusting for bin 
widths). A pseudo-S/N is feature relative, and calculated 
within the maximum boundaries of the feature in question. 

This procedure is repeated 100 times for each noise am- 
plitudcrH For each filter strength and pseudo-S/N we thus 
have a range of pEW errors, from which we obtain the av- 
erage and dispersion. Two sample mappings of these errors 
arc displayed in Figure IB. II (for these maps we have used 
absolute errors). It is seen that for any pseudo-S/N it is 
possible to define filter strengths yielding small errors (dark 
shades in figure) , but the optimal filter strength varies with 
pseudo-S/N. 

These maps are used to find the correct filter for a given 
feature and pseudo-S/N. Separate maps are created for each 
feature, where broad features typically demand stronger fil- 
tering. Furthermore, the dispersion of pEW-values in the 
optimal bin can be used to approximate the systematic er- 
ror of doing pEW measurements on noisy spcctr43- 

Note that it is the shape of the feature that determines 
correct filtering, and that this evolves with epoch. To cor- 
rectly account for this, the above proc edure was repeate d 
for each epoch of the Hsiao templates dHsiao et al.ll2007l ). 
The templates were interpolated to 2.5 A bins in all simu- 
lations. 



Repeated tests were run to verify that results were not sen- 
sitive to the number of iterations. 

12 This systematic error would only include pure noise effects 
and not e.g. effects like host galaxy contamination. 



Simulation results are written to a table. These provide, 
for every feature and lightcurve epoch, the best filter- width 
to use to minimise the risk for noise bias. Since only the 
pseudo-S/N is used, we do not require error spectra. 

The application to real data can be summarised as: 

1. A minimal boxcar is applied, through which the pseudo- 
S/N is determined. 

2. By comparing Monte-Carlo runs for the Hsiao template 
of the same epoch and feature, the optimal boxcar width 
is determined. 

3. The average MC error and dispersion around the ref- 
erence values are taken as systematic errors from the 
simulation. 



B.3. Velocity noise errors 

For the well-defined Type la SN minima studied here, min- 
imum positions are stable relative to noise as long as suffi- 
ciently wide bins are used. A constant bin width in velocity 
space can thus be used. However, determinations of minima 
will still be affected by noise to the degree that on average 
noisy data will have larger dispersion. Both these effects, 
that no bias occurs and the increased dispersion, were stud- 
ied using MC simulations of the Hsiao templates using the 
same approach as for pEWs. Rando m noise is added to the 
Hsiao templates ()Hsiao et al.ll2007l) and the velocities are 
calculated after binning. 

For every template epoch and feature, both bias and 
dispersion are obtained as functions of pseudo-S/N. For ve- 
locities 2, 3, 5, 6 and 7 (and reasonable epoch intervals), 
these results are consistent with no bias and a gradual in- 
crease in dispersion with noise. 

For each spectrum studied (in both the reference and 
NTT/NOT set), epoch and pseudo-S/N values were used 
to locate the corresponding MC dispersion, which was used 
as systematic velocity error. 

For features with more complicated minima (feature 4) 
or possible additional high velocity absorption features (fea- 
ture 1), simply determining the minima will not be enough. 
These features demand either stringent minima criteria or 
function fitting for optimal study. Automatic minima mea- 
surements will show a large scatter. 



Appendix C: Data tables 



28 



Nordin, Ostman et al.: Spectral properties of Type la supernovae up to z ~ 0.3 
Table C.l. Supernova spectra. 



SN 



epochs (days) 



Spec source 



(1992f 



LC source 



SN1983g 6, 7 

SN1986g -4, -4, -4, -3, -3, -2, -1, 0, 
1, 21 

SN1989b 0, 6, 11, 21, 22 
SN1990n 2 

SN1991bg 1, 1, 2, 2, 16, 18, 20, 29 

SN1991m 3, 28 
SN1991s 14 

SN1991t -12, -11, -10, -9, -8, -7, -6 
SN1992g 15 

SN1994d -11, -11, -10, -9, -8, -5, -4, 
-2, 2, 2, 3, 4, 4, 5, 7, 10, 
10, 11, 11, 12, 13, 15, 17, 

19, 24, 26 
SN1994q 10 
SN1994s 22 

SN1996x -4, -2, -1, 0, 1, 3, 7, 12, 22, 
24 

SN1997br -9, -8, -7, -6, -4, 8, 24 
SN1997cn 3, 28 

SN1997do-12, -11, -8, -7, 8, 10, 11, 

12, 14, 15, 20, 21 
SN1997dt -11, -10, -9, -8, -5, 0, 2 
SN1998ab-8, 17, 18, 19, 20, 21, 22 
SN1998aq-9, -8, -3, 0, 0, 1, 1, 2, 2, 3, 
3, 4, 4, 5, 5, 6, 6, 7, 18, 19, 

20, 21, 23, 24, 30 
SN1998bu-4, -4, -3, -3, -2, -2, 0, 8, 8, 

9, 9, 10, 10, 10, 11, 11, 12, 
12, 13, 13, 27, 27, 28, 28, 
29, 29, 30 
SN1998dh-9, -9, -8, -6, -4, -1 
SN1998dm4, 5, 7, 10, 12, 15, 17, 24 
SN1998eg -1, 4, 5, 17, 19, 23 
SN1998v 0, 1, 2, 11, 12, 14 
SN1999aa-ll, -11, -10, -9, -8, -7, -7, 
-7, -6, -5, -4, -3, -3, -3, -2, 
-1, -1, 0, 5, 5, 6, 6, 14, 14, 
14, 15, 16, 17, 19, 19, 25, 
25, 27, 28, 28, 28, 29, 30 
SN1999ac -15, -15, -11, -9, -9, -3, 0, 
0, 2, 2, 2, 7, 
16, 16, 24, 28, 28 

SN1999af -5, 1, 15, 17, 17, 25 
SN1999ao5, 7, 9, 12, 17 
SN1999ar 5 

SN1999aull, 15, 18, 21 
SN1999av 2, 5, 9, 30 
SN1999aw3, 5, 9, 12, 15, 23, 30 
SN1999be 14, 19, 26 
SN1999bi 5, 11, 12, 26 
SN1999bk4, 6, 8 
SN1999bm3, 5, 24 
SN1999bn2, 12, 19, 24 
SN1999bp-2, 0, 1, 6, 16, 21 
SN1999bq3, 3, 16, 20, 24 
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Table C.l. continued. 



SN 



epochs (days) 



Spec source 



LC source 



'Matheson ct al. (2008): Arseniievic ct al. ( 2008) 



SN1999by-5, -5, -4, -4, -3, -3, -2, -2, 

1, 2, 3, 3, 3, 4, 4, 5, 5, 6, 
6, 6, 7, 7, 8, 9, 10, 11, 16, 

24, 25, 28, 29 
SN1999cc -4, -2, -1, 1, 18, 23, 25 
SN1999ee -11, -9, -4, -2, 0, 5, 7, 9, 14, 

17, 20, 25, 30 
SN1999ej -1, 2, 4, 8, 11 
SN1999gd2, 9, 27 
SN1999gp-5, -2, 0, 3, 5, 7, 22 
SN2000cf 3, 4, 14, 16, 24, 25 
SN2000cn -10, -9, -8, 8, 10, 12, 21, 26, 

27 

SN2000cx -4, -3, -2, -1, 0, 1, 5, 6, 7, 9, 
11, 14, 19, 22, 24, 26, 28, 
30 

SN2000dk-5, -4, 1, 4, 9 
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Table C.2. NTT/NOT spectra. 



ID 


IAU 


SPID 


Epochs (days) 


12781 


2006er 


680 


10.9 


12843 


2006fa 


727 


10.2 


12853 


2006ey 


685 


10.3 


12856 


2006A 


695 


-3.2 


12860 


2006fc 


688 


-1.9 


12898 


2006fw 


712 


-6.6 


12930 


2006ex 


687 


10.1 


12950 


2006fy 


700 


-4.4 


13025 


2006fx 


761 


3.4 


13044 


2006fm 


724, 1062 


-8.2, 20.2 


13070 


2006fu 


736 


6.9 


13072 


200611 


723 


0.0 


13135 


2006fz 


739, 998 


-7.7, 17.6 


13796 


2006hl 


1058, 1058 


12.7, 12.7 


13894 


2006jh 


1039 


9.2 


14157 


2006kj 


1040 


9.4 


14437 


2006hy 


1061 


14.0 


14846 


2006jn 


1014 


-1.7 


14871 


2006jq 


1008 


-4.2 


14979 


2006jr 


1009 


-2.1 


14984 


2006js 


1027 


-1.2 


15129 


2006kq 


1015 


1.8 


15132 


2006jt 


1012 


-2.4 


15161 


2006jw 


1010 


-1.0 


15171 


2006kb 


1045, 1045 


-5.7, -5.7 


15203 


2006jy 


1026 


-2.4 


15222 


2006jz 


1004 


-5.8 


15259 


2006kc 


1051 


-1.9 


16021 


2006nc 


1355 


11.3 


16069 


2006nd 


1467 


11.5 


16165 


2006nw 


1326 


2.6 


16215 


2006ne 


1456 


4.3 


16287 


2006np 


1449, 1449, 1569, 


2.4, 2.4, 3.3, 3.3, 






1569, 1650 


19.5 


16352 


2006pk 


1478 


4.1 


16473 


2006pl 


1520 


1.3 


16637 




1514 


-0.9 


17332 


2007jk 


1899 


3.3 


17366 


2007hz 


1782 


8.9 


17389 


2007ih 


1811 


7.0 


17435 


2007ka 


1902, 1902 


2.7, 2.7 


17497 


2007jt 


1837 


-2.4 


17552 


2007jl 


1789 


3.5 


17745 


2007ju 


2161 


15.1 


17784 


2007jg 


1842 


-5.5 


17790 


2007jx 


1887 


1.0 


17811 


2007ix 


1816, 1816 


4.6, 4.6 


17825 


2007jc 


1819 


-4.9 


17875 


2007jz 


1817 


0.3 


17880 


2007jd 


1843, 1957 


-1.9, 1.2 


17886 


2007jh 


1844 


-4.5 


18325 


2007mv 


2277 


8.6 


18466 


20071m 


2270 


4.5 


18768 


20071h 


2135 


6.7 


18787 


2007mf 


2150 


0.2 


18804 


2007me 


2148 


-5.0 


19023 


20071s 


2236 


-1.7 


19101 


2007ml 


2268, 2268 


-6.0, -6.0 


19149 


2007m 


2275, 2275 


-7.1, -7.1 


19155 


2007mn 


2607 


18.7 
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Table C.2. continued. 



ID 


IAU 
SPID 


Epochs (days) 






19282 


2007mk 


2280 


-8.2 




19341 


2007nf 


2298 


-2.4 




19353 


2007nj 


2281 


-7.3 




19381 


2007nk 


2283, 2283 


-3.5, 


-3.5 


19899 


2007pu 


2550 


1.2 




19913 


2007qf 


2585 


9.6 




19953 


2007pf 


2602 


4.2 




19968 


2007ol 


2549 


5.2 




20039 


2007qh 


2584 


7.6 




20040 


2007rf 


2612 


6.4 




20142 


2007qg 


2586 


4.7 




20144 


2007ql 


2541 


1.1 




20345 


2007qp 


2567, 2567 


-0.7, 


-0.7 


20364 


2007qo 


2581 


-1.3 




20430 


2007qj 


2543 


1.4 




20625 


2007px 


2551, 2604 


-5.4, 


-3.6 


21006 


2007qs 


2566 


1.7 




21033 


2007qy 


2565 


-3.3 




21034 


2007qa 


2719 


13.3 




21042 


2007qz 


2564 


-6.4 




21422 


2007rq 


2599 


-3.6 




21502 


2007ra 


2574, 2575 


-8.6, 


-7.7 
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